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1  Executive  Summary 

Our  work  on  this  program  began  in  2004  with  a  long-shot:  We  would  attempt  to  construct 
a  new  £rst-of-kind  source  of  biphotons  and  conditional  single  photons  that  have  a  linewidth 
on  the  scale  of  atomic  transitions,  that  is  about  several  MHz.  Such  photons  are  about 
a  microsecond  long.  A  single  photon  therefore  has  a  physical  length  equal  to  that  of  a 
football  held.  Our  proposed  technique  was  based  on  the  use  of  electromagnetically  induced 
transparency  and  slow  light,  where  basically  the  temporal  length  of  the  biphoton  is  equal 
to  the  physical  length  of  our  atomic  sample  (about  two  cm)  divided  by  the  optical  group 
velocity.  This  idea  worked  almost  immediately  and  the  publication  “Generation  of  Paired 
Photons  with  Controllable  Waveforms,”  appeared  in  Physical  Review  Letters  in  May  of  2005. 

As  time  progressed,  the  objectives  of  this  program,  with  joint  support  from  AFOSR  and 
ARO,  was  aimed  at  both  improving  the  performance  of  this  new  light  source,  and  also  at 
demonstrating  novel  effects  that  are  only  possible  with  long  biphotons.  Work  on  a  new 
quantum  effect  termed  as  nonlocal  modulation  and  requiring  the  use  of  short  biphotons  was 
included  in  the  program. 

We  begin  by  summarizing  a  key  property  of  time-energy  entangled  biphotons.  This  is:  if 
an  observer  at  point  A  chooses  to  measure  the  frequency  of  an  arriving  photon  he  will  then 
know  to  high  accuracy  the  frequency  of  the  photon  which  will  be  measured  by  an  observer 
at  point  B.  But  instead,  if  the  observer  at  point  A  chooses  to  measure  the  time  of  arrival  of 
a  photon  at  his  location,  he  will  then  know,  again  to  high  accuracy,  the  time  of  arrival  of 
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the  photon  at  point  B.  The  accuracy  of  these  measurements  is  not  limited  by  the  Heisenberg 
uncertainty  principle. 

We  turn  next  to  what  is  meant  by  long  and  short.  Typical  biphotons  as  generated  by 
nonlinear  optical  crystals  have  temporal  lengths  in  the  range  of  between  0.1  ps  and  10  ps. 
Such  photons  are  not  resolvable  by  presently  existing  photo  detectors.  These  detectors, 
measure  whether  or  not  a  photon  is  present,  but  may  not  be  used  to  examine  the  functional 
form  of  the  photonic  wave  packet. 

A  key  accomplishment  of  the  present  program  was  the  generation  of  photons  whose  length 
could  be  continuously  varied  from  50  to  900  ns.  Of  importance,  the  line  width  of  these  pho¬ 
tons  is  less  than  the  natural  linewidth  of  the  rubidium  vapor  that  was  used  to  produce  them. 
This  is  important  because  optical  nonlinearities  when  produced  using  electromagnetically 
induced  transparency  continue  to  increase  in  the  subnatural  linewidth  regime  and  in  the 
absence  of  dephasing  may  be  made  arbitrarily  large. 

We  mention  a  surprise  that  occurred  during  this  work.  In  the  course  of  observing  long 
biphotons  with  a  length  determined  by  the  slow  optical  group  velocities,  we  found  that 
the  photonic  wave  packets  had  a  sharp  leading  edge  spike  on  their  front  edge.  Following 
a  suggestion  by  Dan  Gauthier  we  recognized  that  this  front  edge  spike  is  a  Sommerfeld- 
Brillouin  precursor.  This  observation  is  important  because  it  clarifies,  for  both  slow  light 
and  fast  light,  that  information  will  always  be  transmitted  at  the  speed  of  light  in  vacuum. 

In  January  2008  we  recognized  that  we  had  the  capability  to  modulate  single  photons  for 
the  hrst  time.  To  do  this  we  used  the  Stokes  photon  of  a  biphoton  pair  to  set  the  time  origin 
for  electro-optic  modulation  of  the  wave  function  of  the  anti-stokes  photon.  With  the  time 
origin  determined,  the  modulator  could  arbitrarily  modulate  either  the  amplitude  or  phase 
of  the  anti-stokes  photon.  The  technique  therefore  provides  the  technology  for  studying  the 
response  of  atoms  to  shaped  single-photon  waveforms  on  a  time  scale  comparable  to  the 
natural  linewidth  of  target  atoms. 

The  next  step  in  our  work  in  modulating  biphotons  was  the  development  of  a  method  for 
measuring  their  length  using  slow  detectors.  The  essential  idea  is  that  modulation  in  the  time 
domain  followed  by  slow  integration  constitutes  a  Fourier  transformation.  The  experimental 
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technique  is  to  measure  the  coincidence  count  rate  between  single  photon  counting  modules  as 
a  function  of  an  applied  sinusoidal  frequency.  The  inverse  Fourier  transform  of  the  data  then 
yields  the  biphoton  waveform.  Though  this  experiment  was  a  proof  of  principle  experiment, 
ultimately  it  could  be  used  to  measure  wavepacket  profile  of  biphotons  when  sufficiently  fast 
photo  detectors  are  not  available. 

In  recent  months  we  have  demonstrated  an  important  extension  of  our  work  on  modulating 
single  photons.  This  is  the  application  of  spread  spectrum  techniques  at  the  single  photon 
level.  Spread  spectrum  is  well  known  in  the  communications  industry  as  a  technique  for 
avoiding  interference  and  jamming,  and  at  times  increasing  information  capacity.  Our  work 
is  the  hrst  demonstration  of  this  technique  to  single  photons.  As  described  in  this  report  this 
has  allowed  us  to  send  a  single  photon  with  a  known  temporal  shape  through  an  environment 
of  thousands  of  times  more  intense  background  photons  deliberately  produced  by  a  jamming 
laser.  We  have  thus  shown  that  a  single  photon  as  used  for  to  send  a  secure  cryptographic 
key  may  be  deliberately  hidden,  and  recovered  at  the  receiving  end. 

During  the  course  of  this  work  we  recognized  that  there  should  be  a  new  quantum  effect 
that  we  have  termed  as  nonlocal  modulation.  Assume  that  single  and  idler  photons  pass 
through  sinusoidal  phase  modulators  located  at  different  locations.  These  modulators  are 
driven  at  the  same  modulation  frequency  and  are  connected  by  cable  such  that  their  relative 
phase  may  be  varied.  After  passing  through  the  modulators  the  single  and  idler  photons  are 
dispersed,  for  example  by  a  prism,  and  the  relative  positions  of  the  single  and  idler  photons 
are  correlated.  We  find  and  have  experimentally  demonstrated  :  When  the  modulators  are 
run  with  the  same  phase  the  modulation  depths  add;  when  they  are  run  in  phase  opposition 
the  modulation  depths  subtract.  Two  distant  modulators  with  the  same  modulation  depth 
and  opposite  phase  therefore  have  the  same  frequency  correlation  as  when  both  modulators 
are  absent.  This  effect  is  entirely  quantum  mechanical.  Mathematically  it  results  because, 
quantum  mechanically,  one  adds  probability  amplitudes  before  squaring,  while  classically 
one  squares  before  adding. 
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Figure  1:  Energy  Schematic  for  Biphoton  Generation 

2  First  Experiments 

Figure  1  shows  an  energy  schematic  of  the  process  of  the  biphoton  generation  process.  In 
the  presence  of  two  counter-propagating  cw  pump  and  coupling  lasers  paired  spontaneous 
Stokes  and  anti-Stokes  photons  are  generated  and  propagate  in  opposite  directions.  The 
anti-Stokes  photon  has  a  center  frequency  equal  to  that  of  a  highly-absorbing  transition  in 
the  atom.  The  coupling  laser  creates  transparency  and  causes  the  anti-Stokes  wave  to  travel 
at  a  slow  and  variable  group  velocity.  This  group  velocity  determines  the  width  of  the  paired 
photon  wave  packet. 

Our  Erst  experiments  demonstrating  the  generation  of  100  ns  time  scale  biphotons  were 
done  using  a  standard- design  all  metal  spherical  MOT  with  an  atom  cloud  with  a  diameter 
of  a  about  2  mm.  In  order  to  contrast  with  the  results  of  the  greatly  improved  2D  MOT 
described  below,  Figures  2  and  3  show  these  hrst  results.  The  shape  of  the  wave  packet 
depends  on  the  relation  of  two  characteristic  times  that  are  the  essence  of  this  work.  The 
hrst  is  a  Rabi-like  coherence  time  and  the  second  is  the  group  delay  time.  The  experimental 
result  in  the  regime  where  the  Rabi  time  is  less  than  the  group  delay  time  is  shown  in  Fig.  2. 
One  observes  behavior  like  that  of  a  single  atom  where,  at  t  =  0,  a  Stokes  photon  causes  the 
excitation  of  the  upper  state  and  the  probability  amplitude  oscillates  between  the  ground 
and  upper  state.  When  the  group  delay  is  sufficiently  large  that  the  delay  time  is  longer 
than  the  Rabi  time,  then  the  wave  packet  approximates  the  correlation  function  of  an  ideal 
parametric  emitter.  Here,  the  width  of  the  wave  packet  is  approximately  equal  to  the  delay 
time.  Figure  3  shows  experiment  and  theory  as  the  parameters  approach,  but  are  not  quite 
in  the  ideal  gronp  delay  regime. 
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Figure  2:  Generation  in  the  oscillatory  regime 
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Figure  3:  Generation  nearing  the  group  delay  regime 

In  a  second  experiment  still  using  the  spherical  trapped  atom  geometry,  in  May  of  2005 
we  reported  an  improvement  for  generating  narrow  band  time-energy  entangled  photons. 
Here,  working  in  a  standing  wave  right  angle  geometry,  we  use  a  single  driving  laser  with 
two  AOM’s  to  parametrically  pump  the  atoms,  and  also  to  trap  and  cool  the  Rb 

atoms.  This  technique  results  in  a  modification  and  narrowing  of  the  Glauber  intensity 
correlation  function,  and  at  the  same  time  the  elimination  of  Rayleigh  scattering.  An  energy 
level  diagram,  and  the  geometry  and  timing  for  the  experiment  are  shown  in  Fig.  4.  We 
replace  the  three  driving  lasers  of  the  earlier  experiment  with  a  single  (Ti: Sapphire)  laser. 
The  Ti;Sapphire  laser  frequency  is  down  shifted  by  123.5  MHz  so  as  to  be  resonant  with 
the  |5S'i/2,R  =  2)  — |5P3/2,R  =  2)  transition.  This  frequency,  denoted  by  Up,  acts  as 


6 


500  >is  4.5  ms 


(d) 


Trapping  beam 


Figure  4:  Energy  schematic  for  biphoton  generation  using  a  single  pumping  laser,  (a)  Energy 
level  diagram,  (b)  Geometry  (c)  Timing:  The  photon  generation  window  of  500  p.s  is  followed 
by  a  trapping  and  cooling  period  of  4.5  ms.  (d)  Pumping:  The  Ti:Sapphire  laser  is  up-shifted 
to  trap  and  cool  atoms  and  down-shifted  to  serve  as  a  pump  for  the  paired  photon  generation 
process. 

the  coupling  laser  of  earlier  experiments  and  creates  transparency  on  the  |5S'i/2,E  =  1)  — 
I5P3/2,  F  =  2)  transition.  This  same  frequency  pumps  the  four  wave  down-conversion  process 
that  generates  the  paired  photons.  To  cool  and  trap  ®^Rb  atoms  the  incident  Ti:  Sapphire 
laser  beam  is  up-shifted  by  123.5  MHz  to  20  MHz  below  the  |5S'i/2,  P  =  2)  — >  I5P3/2,  P  =  3) 
transition  (not  shown).  This  is  done  periodically  so  as  to  create  a  10%  duty  cycle,  with  the 
trapping  process  occurring  for  4.5  ms,  followed  by  an  experimental  window  of  500  /rs.  The 
trapping  magnetic  held  remains  on  during  this  experimental  window. 

Of  importance  we  hnd  that  the  use  of  a  retro-rehected  laser  for  both  pumping  and  coupling 
results  in  an  unexpected  change  in  the  shape  of  the  Glauber  intensity  correlation  function  for 
the  Stokes  and  anti-Stokes  photons.  In  the  earlier  Stanford  work  the  shape  of  the  correlation 
function  was  set  by  the  longer  of  two  characteristic  times.  These  times  are  the  group  delay 
time  between  the  Stokes  and  anti-Stokes  photons,  and  the  inverse  Rabi  frequency  of  the 
coupling  laser.  Here,  because  of  the  much  larger  Rabi-frequency  (200  MHz)  of  the  common 
pump  and  coupling  laser,  the  group  delay  time  is  short  and  the  Rabi-period  dominates. 
We  therefore  expect  a  damped-  periodic  correlation  function  similar  to  that  of  the  earlier 
work  with  a  period  of  about  5  ns.  We  hnd  experimentally  and  verify  theoretically  that  the 
interference  fringes  (in  the  x-direction)  that  are  caused  by  the  counter-propagating  pumping 
beams  result  in  a  reduction  of  the  outer  Rabi  side  lobes  of  the  intensity  correlation  function. 
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Figure  5:  Coincidence  count  rate  as  a  function  of  the  delay  between  detected  photons,  (a) 
Pump  power  =  160  mW;  flco  =  8O713  and  flpo  =  I4O713.  (b)  Pump  power  =  40  mW, 
flco  =  4O713  and  flpo  =  7O713.  The  pump  is  tuned  to  resonance  with  the  |2)  — >  |3)  transition. 
Data  (o)  was  collected  over  200  s.  Theoretical  curves  (solid  lines)  are  scaled  vertically  by  a 
common  factor  and  are  averaged  over  a  1  ns  bin. 


These  results  are  shown  in  the  Fig.  5  where  we  plot  coincidence  counts  versus  delay  time 
between  the  Stokes  and  anti-Stokes  photons.  Here,  the  pump  laser  is  tuned  to  resonance  of 
the  1 2)  — >  1 3)  transition.  Since  the  Stokes  and  anti-Stokes  photons  are  present  in  both  hbers, 
the  correlation  function  is  symmetric. 


3  Subnatural  Linewidth  Biphoton  Generation  with  a 
2-D  MOT 

A  major  breakout  is  this  program  occurred  with  the  successful  operation  of  a  two-dimensional 
MOT  which  allowed  much  higher  optical  depth  and  performance  in  all  respects.  Of  particu¬ 
lar  importance,  it  allowed,  for  the  hrst  time,  the  generation  of  photons  having  a  sub-natural 
linewidth.  Figure  6  shows  both  experimental  conhguration  (geometry)  and  mechanism  of 
parametric  paired-photon  generation.  In  contrast  to  previous  (spherical)  MOT,  the  new 
2-D  ®^Rb  MOT  has  a  cigar  shaped  atom  cloud(~1.7  cm  long  and  an  aspect  ratio  of  25) 
and  consequently  a  large  optical  depth  in  the  longitudinal  direction;  moreover,  its  cylindri¬ 
cal  quadrupole  trapping  held  results  in  minimal  longitudinal  magnetic  held  gradient  and 
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Figure  6;  Upper:  Biphoton  generation  in  a  double-A  system,  (a)  Experimental  configura¬ 
tion.  Fi  and  F2  are  narrow-band  optical  frequency  filters,  (b)  ®®Rb  energy  level  diagram. 
In  the  presence  of  counter-propagating  pump  (up)  and  coupling  ((Uc)  beams,  Stokes  {ujs  )and 
anti-Stokes,  (cUa^)  photons  are  generated  into  opposing  single-mode  fibers.  Lower:  2-D  MOT 
apparatus.  The  vacuum  cell  is  6  cm  size  ceramic  structured  octagon.  The  cell  is  located  in 
the  middle  of  water  cooled  trapping  coil(racetrack-shaped  cage). 


hence  greatly  reduces  the  inhomogeneous  Zeeman  broadening  of  the  m-states  of  the  5S'i/2 
level.  The  experimental  cycle  comprises  4.5  ms  of  trapping  time  and  0.5  ms  paired  photon 
generation  window.  At  the  end  of  the  trapping  cycle,  the  rubidium  cloud  is  prepared  in 
5S'i/2  level  by  turning  off  the  repumping  laser  0.3  ms  before  turning  off  the  trapping  laser; 
counter-propagating,  circularly  polarized,  cw  pump  (cUp)  and  coupling  (cUc)  lasers  are  subse¬ 
quently  turned  on  and  phase-matched,  paired  Stokes  (cUs)  and  anti-Stokes  (uas)  photons  are 
spontaneously  generated  and  propagate  in  opposite  directions  as  shown  in  the  figure. 

The  optical  depth  of  the  2-D  ^®Rb  MOT  can  be  varied  up  to  62,  which  gives  ns  enough 
parameter  space  to  verify  the  relation  between  the  optical  group  delay  and  the  length  of 
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Figure  7:  Biphoton  wave  packet  data  for  three  slow  group  delay  cases.  The  upper  row  plots 
are  anti-Stokes  FIT  scan  data(o)  and  FIT  ht(blue  curve),  the  lower  row  plots  are  paired- 
photon  coincidence  count  data(-l-)  and  predicted  wave  packet  shape  with  time  bin  width  of  1 
ns  for  800  s(lower  row).  Propagation  delay  of  anti-Stokes  pulse(red  traces  in  inserts  of  lower 
row  plots)  are  also  presented  for  three  slow  group  delay ((r^)  cases.  Fxperimental  parameters 
are:  Left:  {jg  ~  50ns)  OD=7,  flc  =  4.2O713,  hip  =  I.I6713,  and  Ap  =  48.67713.  Middle: 
[Tg  ~  320ns)  OD=53,  flc  =  4.2O713,  hip  =  I.I6713,  and  Ap  =  48. 67713.  Right:  (r^  ~  900ns) 
OD=53,  ric  =  2.35713,  hip  =  I.I6713,  and  Ap  =  48. 67713. 

the  biphoton  waveform.  Figure  7  shows  sets  of  anti-Stokes  FIT  scan  and  paired-photon 
coincidence  counts  for  three  anti-Stokes  FIT  group  delay  cases,  which  is  controlled  by  varying 
the  optical  depth  and  the  coupling  laser  power(rg  ~  {2'^i:i/\VLc\‘^)NaL,  where  NaL  is  the 
optical  depth  and  713  is  the  dephasing  rate  of  level  |3)).  The  trapping  laser  used  for  these 
experimental  runs  has  a  power  of  160  mW,  a  beam  diameter  of  2  cm,  and  is  red  detuned  by 
20  MHz  from  the  |5S'i/2,F  =  3)  — |5P3/2,F  =  4)  transition.  A  repumping  laser  is  locked 
to  the  |5S'i/2,  F  =  2)  — >  I5P3/2,  P  =  2)  transition,  has  a  power  of  80  mW,  and  overlaps  one 
of  six  trapping  beams.  The  pump  laser  is  circularly  polarized  (cr“),  has  a  1/e^  diameter  of 
1.46  mm,  and  is  blue  detuned  from  the  |1)  — >■  |4)  transition  by  146  MHz,  i.e.  Ap  =  48.67713. 
The  coupling  laser  is  circularly  polarized  (cr+),  has  a  1/e^  beam  diameter  of  1.63  mm  and 
is  on  resonance  with  the  |2)  — >  |3)  transition.  The  counter-propagating  pump  and  coupling 
beams  are  collinear  and  set  at  a  2  degree  angle  from  the  longitudinal  axis  of  the  MOT.  The 
Stokes  (a~)  and  anti-Stokes  ((T+)  photons  are  coupled  into  opposing  single  mode  hbers  after 
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passage  through  A/4  wave  plates  and  polarization  beam  splitters  (PBS).  The  Stokes  and 
anti-Stokes  hber  coupling  efficiency  is  70%  and  the  1/e^  waist  diameter  of  their  foci  is  220 
fim. 


Figure  8:  (a)  Measured  correlation  time  vs  measured  anti-Stokes  group  delay.  The  solid  line 
is  a  linear  least  squares  £t.  (b)  Paired  counts  in  a  1  ns  bin  in  800  s  as  a  function  of  the 
optical  depth. 

As  shown  in  Fig.  7,  the  temporal  length  of  biphoton  wavepacket  generated  in  high  OD  case 
are  much  longer(by  more  than  one  order)  than  previously  reported  in  our  phase  I  research. 
It  also  means  that  the  generated  biphoton  has  much  narrowed  bandwidth.  The  predicted 
biphoton  packet  waveforms  in  lower  row  plots  of  Fig.  7  are  computed  with  all  parameters 
obtained  from  the  FIT  measurements  and  vertically  scaled  to  fit  the  experimental  data.  The 
calculated  biphoton  linewidths  are  9.66,  2.36  and  0.75  MHz  respectively.  These  linewidths 
are  comparable  to  the  measured  FIT  bandwidths  and  in  the  latter  two  cases  are  less  than  the 
6  MHz  natural  linewidth  of  Rb  D  line.  The  biphoton  generation  now  is  in  the  linear  group 
delay  regime  where  Tg  >  Tr  and  the  correlation  width  directly  follows  group  delay  time.  This 
is  shown  in  Fig.  8(a). 

Having  taken  into  account  the  hlter  and  etalon  transmissions,  the  fiber  to  hber  coupling 
efficiency,  the  detector  quantum  efficiencies  and  the  duty  cycle.  For  the  conditions  of  Fig.  7, 
we  observe  a  total  of  3213,  31674,  and  22000  paired  counts  in  800  seconds,  which  correspond 
to  generation  rates  of  1275,  12569,  and  8730  pair/s,  respectively.  Higher  generation  rates 
can  be  achieved  by  increasing  the  pump  laser  power.  With  Hp  =  6.88713,  the  paired  photon 
generation  rates  are  4.0  x  10“^  and  9.0  x  10^  pair/s  at  OD=17  and  30,  respectively.  Figure  8(b) 
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shows  that  the  number  of  paired  counts  varies  linearly  with  the  optical  depth.  Though 
the  generation  rate  per  spectral  bandwidth  varies  as  the  square  of  the  optical  depth,  the 
bandwidth  reduces  linearly  with  this  depth  leading  to  the  linear  dependence.  Experimentally, 
under  optimum  conditions,  we  observe  74%  of  the  Stokes  photons  to  be  paired.  We  also 
observe  that  all  of  the  correlation  data  violate  the  Cauchy- Schwarz  inequality  by  as  much  as 
a  factor  of  11600. 

It  is  noticeable  that  there  is  a  sharp  peak  at  the  leading  edge  of  the  correlation  data 
generated  for  high  optical  depth  case  as  shown  in  lower  plots  in  middle  and  right  column  of 
Fig.  7.  This  feature  is  Sommerfeld-Brillouin  precursor  type  as  in  the  case  for  propagating 
classical  wave  packets.  The  physical  picture  is  the  following:  the  detectors  register  biphoton 
coincidence  counts  versus  the  time  r  =  tas  —  tg]  the  earliest  portion  of  the  biphoton  wave 
packet  comes  from  the  high  frequency  portion  of  the  spectrum,  which  is  not  in  the  range  of 
large  group  delay,  and  is  thought  of  as  the  Sommerfeld  precursor.  At  slightly  later  times  the 
low  frequency  Brillouin  components  arrive  at  the  detector  and  beat  with  the  simultaneously 
arriving  high  frequency  components.  Though  precursors  are  now  understood  in  the  optical 
region  and  have  even  been  observed  long  ago  with  correlated  gamma-ray  photons,  our  work 
reports  the  hrst  observation  of  precursors  as  measured  by  single  photon  correlation. 

For  more  details  please  see:  1)  Shengwang  Du,  Pavel  Kolchin,  Chinmay  Belthangady, 
G.  Y.  Yin,  and  S.  E.  Harris,  ’’Subnatural  Linewidth  Biphotons  with  Controllable  Temporal 
Length,”  Phys.  Rev.  Lett.  100,  183603  (May,  2008);  2)  Shengwang  Du,  Pavel  Kolchin, 
Chinmay  Belthangady,  G.  Y.  Yin,  and  S.  E.  Harris,  ’’Observation  of  Optical  Precursors  at 
the  Biphoton  Level,”  Optics  Letters  33,  2149  (September,  2008). 
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4  Theory  of  EIT-  based  Paired  Photon  Generation 


Our  theoretical  model  is  based  on  paired  photon  generation  in  the  double-lambda  atomic 
system.  We  hrst  reexamine  the  conditions  required  for  the  system  to  operate  in  the  group 
delay  regime.  We  hnd  that  the  optical  depth  of  the  atomic  system  has  to  be  sufficiently  high 
in  order  to  avoid  the  hltering  of  the  generation  bandwidth  of  paired  photons  by  the  trans¬ 
parency  window.  Second,  we  have  extended  the  theoretical  treatment  beyond  the  ground 
state  approximation.  This  allows  us  to  properly  include  and  analyze  the  effect  of  Langevin 
noise  ffuctuation  on  the  atomic  system  and  solve  the  problem  of  its  return  to  the  ground  state 
after  the  emission  of  Stokes  and  anti-Stokes  photon  pairs.  We  also  addressed  the  important 
questions:  1)  ‘What  are  the  effects  of  the  optical  thickness  of  the  atomic  sample  on  paired 
and  single  photon  generation?”  2)  “Does  every  Stokes  photon  have  its  paired  anti-Stokes 
photon?” 

With  low  parametric  gain  and  high  optical  depth  we  show  that  the  system  can  produce 
highly  correlated  photon  pairs.  The  shape  of  the  intensity  correlation  function  and  the 
emission  bandwidth  depend  on  the  coupling  laser  Rabi  frequency  and  the  optical  depth  of  the 
atomic  sample.  Compared  to  SPDC,  paired  photon  generation  in  the  double-lambda  atomic 
system  is  affected  by  Raman  gain  in  the  Stokes  channel  and  EIT  in  the  anti-Stokes  channel. 
EIT,  through  the  absorption  at  the  poles,  cuts  the  emission  bandwidth.  In  order  to  enter 
a  regime  where  the  EIT  window  is  sufficiently  large  and  therefore  the  emission  bandwidth 
is  controlled  to  a  large  extent  by  the  phase-matching  process  in  the  presence  of  large  group 
delay,  the  optical  depth  of  the  atomic  sample  has  to  be  larger  than  10.  High  optical  depth 
substantially  reduces  the  influence  of  Langevin  noise  fluctuations  and  Raman  scattering  on 
paired  photon  generation  so  that  the  Stokes  and  anti-Stokes  photons  are  generated  mostly 
in  pairs. 

For  more  details  please  see:  P.  Kolchin,  “Electromagnetically-Induced- Transparency-Based 
Paired  Photon  Generation,”  Phys.  Rev.  A  75,  033814  (March  2007).  Kolchin’s  dissertation 
chapter  4  “Theory  of  EIT  based  Paired  Photon  Generation”  is  included  in  this  report  as 
Appendix  C. 
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5  Electro-Optic  Modulation  of  Single  Photons 


A  highlight  of  our  work  on  this  DARPA  program  has  been  the  demonstration  of  conditional 
shaping  of  single  photon  waveforms  using  electro-optic  modulators.  We  use  the  Stokes  photon 
of  a  biphoton  pair  to  set  the  time  origin  for  electro-optic  modulation  of  the  wave  function  of 
the  anti-Stokes  photon.  This  technique  allows  arbitrary  control  of  both  phase  and  amplitude 
of  single  photons.  A  key  requirement  for  this  modulation  scheme  is  that  the  temporal  length 
of  biphoton  wavepacket,  which  is  represented  by  correlation  time  in  experiment,  is  much 
longer  than  modulation  response  time  designed  for  experiment.  The  subnatural  linewidth 
biphoton  source  described  in  the  previous  section  producing  biphotons  with  correlation  times 
adjustable  in  the  50-900  ns  range  is  ideally  suited  for  this  application.  Our  single  photon 
counting  modules(about  350  ps  of  timing  resolution)  and  data  system  has  fast  enough  tem¬ 
poral  resolution  to  observe  AM  modulation  in  this  work. 


Figure  9:  Schematic  of  paired  photon  generation  and  conditional  modulation.  A  Stokes 
photon  detected  by  an  SPCM  sets  the  time  origin  for  shaping  the  anti-Stokes  photon  with 
an  electro-optic  modulator.  To  within  the  accuracy  of  the  SPCM,  this  allows  shaping  of 
both  the  amplitude  and  phase  of  the  anti-Stokes  photon. 


Figure  9  shows  the  schematic  of  the  experiment.  We  use  counterpropagating  cw  pump 
and  coupling  lasers  to  generate  time-energy  entangled  pairs  of  Stokes  and  anti-Stokes  photons 
which  propagate  in  opposite  directions  and  are  collected  into  single  mode  hbers  as  described 
in  section  A.  The  detection  of  the  Stokes  photon  at  D1  sets  the  time  origin  for  bring  the 
function  generator  that  drives  the  electro-optic  modulator  which  in  turn  modulates  the  anti- 
Stokes  photon.  Verihcation  of  the  single  photon  nature  of  the  modulated  anti-Stokes  photon 
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(o') 

is  done  using  a  50-50  beam  splitter  and  detectors  D2  and  D3.  Their  coincidence  data 

is  expected  to  show  a  dip  at  r  =  0  proving  nonclassical  nature  of  modulated  anti-Stokes 

signal. 

The  electro-optic  amplitude  modulator  consists  of  phase  modulators  in  both  arms  of  a 
Mach-Zehnder  (MZ)  interferometer.  The  degree  of  phase  control  in  both  arms  depends  on 
the  type  of  the  modulator.  We  use  a  z-cut  modulator  that  requires  14  =  1-75  volts  to 
cause  the  vr  phase  shift  required  to  go  from  minimum  to  maximum  transmission  and  can 
be  operated  at  a  maximum  frequency  of  10  GHz.  One  port  of  the  output  beam  splitter 
of  the  MZ  interferometer  is  terminated  so  that  the  portion  of  the  photon  wave  function 
that  is  not  transmitted  is  lost.  In  general,  if  a  Stokes  photon  is  detected  at  time  ti,  and 
the  modulator  is  activated  conditioned  on  this  detection  then,  in  the  Heisenberg  picture, 
the  anti-Stokes  operator  at  the  output  of  the  modulator  is  related  to  the  input  operator 
by  ttoutih)  =  f  5'(442)®m(^2)'^^2-  If  there  are  no  dispersive  elements,  then  to  within  an 
unimportant  phase  factor  we  may  write  aout(t2)  =  ■  The  correlation  function  in 

the  presence  of  the  modulator  is  related  to  that  in  the  absence  of  the  modulator  by 

Gg)(r)  =  |m(r)|W^>(r)  (5,1) 

With  the  biphoton  wavefunction  given  by  -|-  r)  the  modulated  (conditional)  single 

photon  wavefunction  is  m(r)T(fi,fi  -|-  r).  We  adjust  the  bias  voltage  at  the  input  of  mod¬ 
ulator  so  that  the  output  of  the  modulator  m(r)  is  related  to  the  input  voltage  V{t)  by 
m(r)  =  sin  [(/)(r)]  exp  [ia(/)(r)],  where  0(r)  =  7rl/(r)/(214),  and  a  is  a  phase  modulation 
parameter.  For  a  z-cut  amplitude  modulator  as  used  here  a  =  0.75,  but  may  be  eliminated 
by  using  an  x-cut  modulator. 

The  principal  experimental  results  of  this  work  are  shown  in  Figure  10.  In  part  (a),  the 
modulation  signal,  shown  as  an  inset,  is  a  set  of  two  square  pulses.  Of  importance,  there  is 
no  vertical  scaling  between  the  modulated  and  non- modulated  waveforms.  In  Fig.  10(b), 
we  show  photons  modulated  with  two  different  waveforms.  In  the  hrst  case  the  modulator 
is  driven  with  a  gaussian  pulse.  In  the  second  case  we  design  the  function  generator  output 
so  as  to  compensate  for  the  nonlinear  distortion  of  the  modulator  in  such  a  way  that  the 
output  of  the  modulator  is  an  exact  rising  exponential. 
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Figure  10:  D1-D2  coincidence  counts  in  a  1  ns  bin  as  a  function  of  the  delay  between 
Stokes  and  anti-Stokes  photons,  (a)  Modulated  (a)  and  unmodulated  (□)  waveforms,  (b) 
Waveforms  with  Gaussian  (□)  and  rising  exponential  (A)  shapes.  The  experimental  data  (□, 
a)  were  collected  over  2000  s.  The  solid  curves  for  cases  (a)  and  (b)  are  plotted  from  theory. 
The  inset  in  part  (a)  is  the  scope  trace  of  the  output  voltage  of  the  function  generator. 

We  dehne  the  retrieval  efficiency,  Er  of  a  paired  photon  source  as  the  probability  to  gener¬ 
ate  a  single  anti-Stokes  photon  on  the  condition  that  its  paired  Stokes  photon  is  detected.  For 
the  non-modulated  photon  we  measure  Sr  =  3.5%.  When  losses  at  the  beamsplitter,  modu¬ 
lator,  hlters,  hber  to  hber  coupling  and  detector  efficiency  are  backed  out,  this  corresponds  to 
a  retrieval  efficiency  of  55%.  For  the  modulated  photons  the  measured  retrieval  efficiencies  of 
the  two  square  pulses,  the  rising  exponential  and  gaussian  waveforms  are  £r  =  1.3%,  0.61% 
and  0.9%  respectively.  With  losses  backed  out  these  efficiencies  are  21%,  9.4%  and  11.2% 
respectively. 

Since  single  photons  incident  on  a  beamsplitter  must  go  into  one  output  port  or  the  other, 
in  the  ideal  case  where  there  are  no  two-photon  events  and  there  is  no  excess  scattered  light. 
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Figure  11:  (color  online)  Conditional  three-fold  correlation  function  as  a  function 

of  the  Stokes  rate  for  unmodulated  (□)  and  modulated  (A)  single-photon  generation.  The 

f  2') 

dashed  curve  shows  the  theoretical  limit  for  fi'cond(O)  absence  of  excess  light  scattering 

(see  text). 

we  would  expect  no  three-fold  coincidences  at  the  detectors.  A  measure  of  the  quality  of 
heralded  single  photons  that  quantihes  suppression  of  two  photon  events  is  given  by  the 
conditional  correlation  function: 

£LM  =  (5.2) 

2Vi2iVi3 

Here  Ni  is  the  number  of  the  Stokes  counts  at  Hi,  A^i2,  and  Ni^  are  the  number  of  two- fold 
coincidence  counts  within  a  time  window  Tc  at  detectors  Hi,  H2  and  Hi,  H3;  and  A^i23  is 
the  number  of  three-fold  coincidence  counts  within  this  same  time  window. 

In  Figure  11  triangles  and  squares  show  measured  versus  Stokes  rate  with  and 

without  modulation.  The  modulation  is  done  with  the  same  signal  as  in  Fig.  10(a).  We 
set  Tc  equal  to  the  nominal  length  of  the  unmodulated  biphoton  (285  ns).  At  a  Stokes  rate 
of  2.2  X  10“^  sec“^  which  corresponds  to  Hp  =  0.26  F3,  we  obtain  g^cJndi^)  =  0-2  =*=  0.04  and 
5'cond(0)  =  0.21  ±0.07  for  the  unmodulated  and  modulated  waveforms  respectively.  The  fact 
that  the  measured  5'cond(0)  is  less  than  0.5,  (the  limiting  value  for  a  two  photon  Fock  state), 
is  indicative  of  the  near-single  photon  character  of  the  light  source. 
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Because  there  is  a  small  probability  for  the  parametric  down  conversion  process  to  generate 
multiple  pairs  of  biphotons,  even  in  the  absence  of  spurious  light  scattering,  the  conditional 
correlation  function  is  not  zero.  The  dashed  curve  In  Fig.  11  shows  the  theoretical  prediction 
for  the  conditional  correlation  function  that  results  from  such  multiple  scattering  events. 
Because  of  light  scattering  from  both  the  pump  and  coupling  lasers,  the  experimental  curves 
lie  above  this  limiting  value. 

We  perform  two  control  experiments:  In  the  hrst  we  remove  the  30  m  long  optical  hber 
so  as  to  modulate  the  uncorrelated  background  noise  in  the  tail  of  the  correlation  function. 
Here,  we  measure  gcoiidi^)  —  In  fhe  second  experiment  we  apply  modulation  at  random 
times,  using  an  external  10  MHz  digital  signal  as  a  trigger  for  the  function  generator.  As 
expected,  we  observe  a  reduced  rate  of  paired  counts  and  no  change  in  the  shape  of  the 
correlation  function. 

The  method  demonstrated  in  this  work  might  be  used  to  optimally  load  a  single  photon 
into  an  optical  cavity,  or  instead,  to  study  the  transient  response  of  atoms  to  different  single 
photon  waveforms.  In  the  context  of  light-matter  interfaces,  it  may  improve  the  efficiency 
of  storage  and  retrieval  of  single  photons  in  atomic  ensembles.  For  quantum  information 
applications,  both  amplitude  and  phase  modulators  could  be  used  to  allow  full  control  over 
the  single  photon  waveforms.  For  example,  one  could  a  construct  a  single  photon  waveform 
that  is  a  train  of  identical  pulses  with  information  encoded  into  the  relative  phase  differ¬ 
ence  between  consecutive  pulses.  The  importance  of  the  electro-optic  method  is  its  speed 
and  ability  to  modulate  phase  as  well  as  amplitude.  The  technique  provides  the  technol¬ 
ogy  for  studying  the  response  of  atoms  to  shaped  single  photon  waveforms  on  a  time-scale 
comparable  to  the  natural  linewidth. 

For  more  details  please  see:  P.  Kolchin,  Chinmay  Belthangady,  Shengwang  Du,  G.Y.  Yin, 
and  S.E.  Harris,  “Electro-Optic  Modulation  of  Single  Photons,”  Phys.  Rev.  Lett.  101, 
103601  (September,  2008). 


18 


6  Modulation  and  Measurement  of  Time-Energy  En¬ 
tangled  Photons 

In  this  section  of  the  final  report,  we  extend  our  work  on  the  modulation  of  single  photons  to 
the  modulation  of  biphotons.  The  essential  difference  is  that  we  no  longer  measure  and  set 
the  timing  by  the  measurement  of  a  single  photon.  Instead,  biphoton  pairs  arrive  at  random 
at  their  respective  detectors.  A  highlight  of  this  work  is  the  invention  and  first  experimental 
demonstration  of  a  Fourier  technique  that  allows  the  measurement  of  fast  biphotons  using 
slow  detectors.  Figure  12  shows  (a)  the  schematic  of  our  proof-of-principal  experiment  and 
(b)  the  laser/rubidium-atom  interaction  diagram  for  paired  photon  generation.  Parametri¬ 
cally  down-converted  spontaneous  signal  and  idler  photons,  or  as  in  the  experiment  of  this 
work,  Stokes  and  anti-Stokes  photons,  are  incident  on  synchronously  driven  sinusoidal  am¬ 
plitude  modulators.  Without  any  modulation,  in  general,  the  setup  will  directly  yield  the 
Glauber  correlation  function  where  r  is  the  relative  arrival  time  of  the  signal  and 

idler  photons,  with  its  time  resolution  limited  by  the  speed  of  photon  detector  (SPCM)  and 
TDC  system.  As  shown  below,  our  new  method,  by  adding  AM  photon  modulation,  will  give 
temporal  biphoton  correlation  measurement  with  time  resolution  beyond  the  SPCM-TDC 
limit . 

When  amplitude  modulators  mi{t)  and  m2{t)  are  introduced  between  the  down  conversion 
source  and  the  detectors,  the  modulated  correlation  function  is  then  written  as 

+  r)  =  \m2{t  +  T)f  (6.1) 

where  the  subscripts  M  and  0  indicate  cases  of  modulation  on  or  off.  In  real  experiments, 
photons  arrive  at  a  random  time  t  which  is  averaged  out  over  the  data  collection  period  T: 

cfF)  =  M(T)Gf\T), 

><(’■)  =  ;);/  \m.2(t  +  t)\^  dt,  (6.2) 

where  A4(t)  is  the  intensity  correlation  function  of  the  modulators  in  the  signal  and 
idler  channels.  If  both  channels  are  modulated  by  sinusoidal  amplitude  modulators  with 
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frequency  ui  and  a  common  phase  yj,  that  is  by  modulators  rriiit)  =  m2{t)  =  cos(a;t  +  (p), 
then,  irrespective  of  this  phase,  M(t)  =  1/4  +  1/8  cos(2a'r). 

If  the  detectors  are  slow  in  the  sense  that  they  integrate  over  the  length  of  the  bipho¬ 
ton  wave  packet  ,  but  short  as  compared  to  the  inverse  rate  of  biphoton  generation,  the 
measurement  becomes  an  integral 


(6.3) 


We  neglect  the  DC  term  and  normalize  to  obtain  the  Fourier  cosine  transform  pair 


(6,4) 


In  the  measurement  procedure,  F{2uj)  is  the  measured  coincidence  count  rate  between  low 
speed  single  photon  counting  modules  (SPCMs)  as  a  function  of  the  sinusoidal  modulation 
frequency  uj.  The  slow  detection  system  (SPCM-I-TDC)  means  that  hardware  integrates 
signals  over  r,  which  is  the  relative  arrival  time  of  the  signal  and  idler  photons.  We  then 
apply  the  above  inverse  Fourier  cosine  transform  to  yield  the  Glauber  correlation  function 
Gq^\t)  and  therefore  the  square  of  the  absolute  value  of  the  biphoton  wavefunction. 

Our  experiment  makes  use  of  long  biphotons  that  are  produced  using  the  techniques  of 
electromagnetically  induced  transparency  and  slow  light.  The  use  of  long  biphotons  allows 
us  to  compare  the  correlation  function  measured  by  our  Fourier  transform  technique  with 
a  direct  measurement  using  fast  detectors  and  a  TDC.  The  experimental  conhguration  is 
shown  in  Fig.  12(a).  Paired  photons  are  generated  with  cold  Rb  atoms  using  the  method 
of  Balic  et  ah  We  apply  strong  counterpropagating  pump  and  coupling  lasers  (not  shown) 
to  produce  phase  matched  counter-propagating  pairs  of  time-energy  entangled  Stokes  and 
anti-Stokes  photons.  We  use  a  ®^Rb  two-dimensional  magneto-optic  trap  with  an  optical 
depth  which  can  be  varied  between  10  to  60  to  generate  biphotons  with  temporal  lengths 
between  50  and  900  ns.  The  inset  in  Fig.  13(a)  shows  the  biphoton  wavefunction  obtained  at 
an  optical  depth  of  35  as  directly  measured  using  a  TDC.  Two  features  are  of  interest.  First, 
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5Si/2  (F=2) 


Figure  12:  (a)  Schematic  of  experiment.  Long  biphotons  are  generated  in  the  cold-atom  Rb 
cell.  Using  an  optical  hber,  the  anti-Stokes  photon  is  delayed  by  175  ns,  and  the  photons  are 
modulated  by  synchronously-driven  sinusoidal  modulators  before  correlation  detection  with 
a  time-to-digital  converter,  (b)  Energy  level  diagram  for  paired  photon  generation  in  Rb. 

the  width  of  the  wavefunction  is  determined  by  the  slow  group  velocity  of  the  anti-Stokes 
photon  and  varies  linearly  with  the  optical  depth.  Second,  the  distinctive  sharp  feature  at 
the  leading  edge  is  a  Sommerheld-Brillouin  precursor  that  ensures  that  the  earliest  signal 
reaches  a  detector  at  the  speed  of  light  in  vacuum. 

The  generated  Stokes  and  anti-Stokes  photons  are  transmitted  through  10  GHz  electro¬ 
optic  amplitude  modulators  (Eospace  Inc.)  with  a  half-wave  voltage,  K-  of  1.3V.  To  obtain 
a  perfectly  sinusoidal  output,  the  modulators  are  biased  at  maximum  transmission  and  the 
input  voltage  is  varied  linearly  using  a  triangular  waveform  that  varies  between  — IG  and  -|-Ur- 
This  waveform  is  generated  by  a  fast  function  generator  (Tektronics  AFG3252)  with  two 
output  channels  whose  frequencies  and  phases  can  be  varied  independently.  The  modulated 
photons  are  then  sent  to  SPGMs  (Perkin  Elmer  SPGM-AQR-13),  which  are  connected  to  the 
start  and  stop  inputs  of  a  TDG  (Fast-Gomtec  TDG  7886S).  Goincidence  counts  are  binned 
into  histograms  and  plotted  as  a  function  of  the  time  difference  between  the  detection  of  a 
Stokes  and  an  anti-Stokes  photon. 

In  Fig.  13  (color  online)  we  demonstrate  the  modulation  of  a  biphoton  wavefunction.  The 
data  are  recoded  by  binning  coincidence  counts  versus  time  into  1  ns  bins.  In  Fig.  13(a),  the 
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Figure  13:  (color  online)  Modulation  of  the  biphoton  correlation  function,  (a)  rriiit)  =  1  and 
m2{t)  =  cos(ci;t  +  (p).  The  inset  in  (a)  shows  the  correlation  function  with  both  modulators 
open,  (b)  mi(t)  =  m2(t)  =  cos(a;t  +  p).  Here  u  =  27rx  35  xlO®. 

Stokes  modulator  is  turned  off  and  is  biased  at  maximum  transmission;  and  the  anti-Stokes 
modulator  is  driven  at  35  MHz.  In  agreement  with  theory,  the  biphoton  wavefunction  is  not 
modulated.  In  Fig.  13(b),  both  modulators  are  modulated  at  35  MHz  with  the  same,  but  ar¬ 
bitrary,  phase.  The  correlation  function  is  now  modulated  at  twice  the  applied  frequency,  i.e. 
at  70  MHz.  When  driven  by  non-sinusoidal  waveforms,  the  correlation  function  is  modulated 
by  the  cross-correlation  of  the  two  modulating  signals.  In  agreement  with  Eq.  (6.2),  we  have 
verihed  that  when  the  two  modulating  signals  are  square  waves  with  the  same  frequency,  the 
correlation  function  is  modulated  by  a  triangular  function. 

We  next  demonstrate  the  Fourier  measurement  technique.  A  35  m  long  polarization- 
maintaining  hber  is  used  to  delay  the  anti-Stokes  photon  by  175  ns  in  the  experiment.  The 
modulators  are  driven  synchronously  at  frequencies  between  0  and  30  MHz.  To  simulate 
slow  detectors,  we  bin  coincidences  into  1  /is  bins,  and  since  the  temporal  extent  of  the 
waveform  is  less  than  1  ps,  the  entire  biphoton  is  contained  within  the  first  time  bin.  From 
the  Fourier  transform  property  that  delay  in  the  time  domain  corresponds  to  oscillation  in 
the  frequency  domain,  the  plot  of  coincidence  counts  versus  modulation  frequency  shows 
ripples  with  a  frequency  (5.7  MHz)  equal  to  the  inverse  of  the  time  delay.  We  show  two 
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Figure  14:  Fourier  Transform  measurement  technique  at  an  optical  depth  of  15 (Left)  and 
35(Right).  Upper:  Frequency  domain  data,  Lower:  Fourier  Cosine  transforms  (red)  of 
corresponding  frequency  domain  data,  and  real-time  fast  direct  temporal  correlation  data 
(blue).  The  term  modulation  frequency  on  the  x-axis  of  (a)  and  (c)  refers  to  the  applied 
frequency.  The  observed  modulation  frequency  [Eq.  (6.4)]  is  a  factor  of  2  higher. 

data  sets  in  Fig.  14  for  optical  depth  of  15  (left)  and  35  (right),  for  relatively  short  and 
long  biphoton  wave  packets,  respectively.  In  Figs.  14(b)  and  (d),  we  compare  the  biphoton 
wavefunction  as  directly  measured  with  1  ns  bins  (blue  line)  to  the  Fourier  cosine  transform 
of  the  traces  in  Figs.  14(a)  and  (c),  (red  line),  vertically  scaled  to  match  the  peak  value.  We 
hnd  reasonable  agreement  between  the  two  methods.  We  note  that  the  sharp  spike  in  the 
temporal  trace  near  f  =  0  results  from  the  DC  component  in  the  frequency  domain  trace. 

We  have  shown  how  biphotons  may  be  modulated,  and  how  this  modulation  may  be  used 
to  measure  the  magnitude  of  the  biphoton  wavefunction.  Though  we  have  used  long  bipho¬ 
tons  and  low  modulation  frequencies,  we  believe  that  this  technique  should  be  extendable 
to  short  biphotons.  A  commercially  available  telecommunication  modulator  driven  at  a  fre¬ 
quency  of  60  GHz  and  therefore  modulating  at  120  GHz  will  allow  measurement  of  biphotons 
with  a  minimum  length  of  about  8  picoseconds.  This  is  about  a  factor  of  hve  faster  than 
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state-of-the-art  commercial  SPCMs.  If  the  comparison  is  made  on  the  basis  of  state-of-the- 
art  polymer  light  modulators  operating  at  200  GHz  and  therefore  modulating  at  400  GHz, 
then  the  Fourier  technique,  will  allow  measurement  down  to  about  2.5  ps.  This  is  about 
eight  times  faster  than  the  fastest  reported  superconducting  detector.  Looking  further  to 
the  future,  all-optical  light  modulators  have  been  demonstrated  at  frequencies  greater  than 
1  THz;  thereby  in  principle  allowing  measurement  of  photons  on  femtosecond  time  scales. 
We  note  that  an  advantage  of  the  Fourier  technique  as  compared  to  either  sum  frequency 
correlation  or  Hong-Ou-Mandel  interference  is  that  the  signal  and  idler  photons  may  be  cor¬ 
related  at  distant  detectors  and  do  not  need  to  be  brought  together  at  a  summing  crystal  or 
at  a  beam  splitter. 

For  more  details  please  see:  Ghinmay  Belthangady,  Shengwang  Du,  Ghih-Sung  Ghuu,  G.Y. 
Yin,  and  S.E.  Harris,  “Modulation  and  measurement  of  time-energy  entangled  photons”, 
Phys.  Rev.  A,  Rapid  Gommunications,  80,  031803  (September  2009). 


7  Observation  of  Nonlocal  Modulation  with  Entangled 
Photons 

The  idea  of  what  we  term  here  as  nonlocal  modulation  was  published  about  two  years  ago; 
i.e.,  S.  E.  Harris,  “Nonlocal  Modulation  of  Entangled  Photons,”  Phys.  Rev.  A  78,  021807(R) 
(2008).  This  section  of  the  report  describes  the  hrst  experimental  observation  of  this  effect. 

When  the  photons  of  a  time-energy  entangled  pair  are  sent  through  different  channels 
having  arbitrary  dispersions,  the  dispersion  in  one  channel  may  be  negated  by  dispersion  of 
the  opposite  sign  in  the  other  channel.  This  effect  results  from  a  quantum  mechanical  inter¬ 
ference  and  has  no  classical  analog.  This  is  now  termed  as  nonlocal  dispersion  compensation. 
We  report  the  hrst  observation  of  a  time- frequency  analog  to  nonlocal  dispersion  cancella¬ 
tion  and  term  this  effect  as  nonlocal  modulation.  Gonsider  a  simplihed  concept  system  as 
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shown  in  left  portion  of  Fig.  15,  where  a  monochromatic  pump  generates  non-degenerate 
time-energy  entangled  photon  pairs.  The  signal  and  idler  photons  (Channels  1  and  2)  pass 
through  sinusoidal  phase  modulators.  These  modulators  are  driven  at  the  same  modulation 
frequency,  and  their  relative  phase  may  be  varied.  After  passing  through  the  modulators, 
the  signal  and  idler  photons  are  dispersed,  for  example,  by  a  prism,  and  the  relative  po¬ 
sitions  (frequencies)  of  the  signal  and  idler  photons  are  correlated.  When  the  modulation 
frequency  is  small  as  compared  to  the  spectral  bandwidth  of  the  signal  or  idler,  we  find  a 
consequence  of  time-energy  entanglement  that  we  term  as  nonlocal  modulation.  Specifically, 
these  distant  modulators  act  cumulatively  to  determine  the  apparent  modulation  depth.  If 
the  two  identical  modulators  have  opposite  phase,  they  negate  each  other  and  act  as  if  nei¬ 
ther  modulator  were  present.  Conversely,  if  operated  with  the  same  phase,  they  produce  the 
same  correlation  as  does  a  single  modulator  with  twice  the  modulation  depth  acting  on  only 
one  of  the  photons. 


Figure  15:  Nonlocal  modulation.  Signal  and  idler  photons  are  phase  modulated  at  the  same 
frequency  and  with  controllable  phase.  The  signal  and  idler  beams  are  frequency-dispersed, 
and  the  positions  (frequencies)  of  the  detected  photons  on  the  photodetectors  are  correlated. 
Left:  simplified  concept  diagram.  Right:  experimental  configuration. 

If  the  resolution  of  the  frequency  correlator  of  Fig.  15  (left)  is  inhnite,  the  resulting 
correlation  traces  for  phased  and  anti-phased  modulators  may  look  something  like  Fig.  16. 
Since  the  pump  laser  which  generates  the  entangled  photons  is  monochromatic,  they  are 
delta-function  correlated  in  frequency  before  any  modulation  occurs.  The  modulators  act 
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Figure  16:  Example  frequency  correlation  function  with  (a)  phased  and  (b)  anti-phased 
sinusoidal  modulators.  Because  the  modulators  have  the  same  amplitude,  cumulative  mod¬ 
ulation  results  in  either  (a)  doubling  of  modulation  depth  or  (b)  cancellation  of  modulation 
in  the  spectral  correlation. 

on  this  delta  function  correlation  by  generating  sidebands  whose  distribution  depends  on 
the  amplitude  of  the  sinusoidal  driving  signal.  If  the  modulators  are  run  in  phase  with 
each  other,  the  distribution  of  sidebands  resembles  that  of  a  single  modulator  with  twice 
the  driving  amplitude  [Fig.  16(a)].  If  run  with  opposite  phase,  the  modulators  cancel  each 
other,  and  a  single  delta  function  is  present  in  the  correlation,  as  if  neither  modulator  were 
present  [Fig.  16(b)]. 

The  more  detailed  experimental  configuration  is  shown  in  right  part  of  Fig.  15.  We 
pump  a  20  mm  long,  periodically-poled,  magnesium  oxide-doped  stoichiometric  lithium  tan- 
talate  crystal  (PPSLT)  with  0.8  W  from  a  532  nm  cw  laser.  The  nonlinear  crystal  is  phase 
matched  to  produce  32  nm  bandwidth,  degenerate  photon  pairs  at  1064  nm.  All  helds  are 
polarized  along  the  extraordinary  axis  of  the  crystal.  The  generated  photons  are  hltered 
from  the  strong  532  nm  pumping  beam  using  a  four-prism  setup  and  are  then  coupled  into 
a  polarization-maintaining  fused-hber  beam  splitter  (AFW  PFC-64-1-50-L-P-7-1-F)  which 
diverts  the  photons  into  Channels  1  and  2  with  equal  probability.  The  photons  are  passed 
through  identical  sinusoidal  phase  modulators  (EOSPACE)  which  are  driven  at  30  GHz  with 
modulation  depths  of  about  1.5  radians.  To  set  the  modulation  depth,  we  adjust  the  variable 
attenuators  in  the  microwave  driver  circuit  and  verify  the  depth  by  measuring  the  sideband 
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amplitudes  of  a  1064  nm  Nd:YAG  reference  laser.  The  relative  phase  between  the  modu¬ 
lators  is  controlled  using  a  calibrated  phase  trimmer.  The  microwave  circuit  schematic  is 
detailed  in  Fig.  17. 
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Figure  17:  Driver  circuit  for  30  GHz  phase  modulators:  (a)  Microwave  Dynamics  PLO-3072 
30  GHz  source,  (b)  Weinschel  Associates  953K-3dB  variable  attenuator,  (c)  Picosecond  Pulse 
Labs  5350-218  power  divider,  (d)  Weinschel  Associates  953K-10dB  variable  attenuator,  (e) 
Atm  Inc.  P1409-360  phase  trimmer,  (f)  Nextec-RF  NA00435  amplifiers,  (g)  MegaPhase  CA- 
V1K2  K  to  V  coaxial  adapters,  (h)  EOSPACE  PM-DV5-40-PFU-PFU-106-LV-UL  custom 
phase  modulators. 


The  photons  then  pass  through  identical  monochromators,  each  having  a  linear  dispersion 
of  210  GHz/mm  and  a  Gaussian  instrument  response  function  with  a  FWHM  of  8.5  GHz.  The 
monochromators  each  consist  of  a  1200  grooves/mm,  aluminum-coated  grating  (Thorlabs 
GR50-1210)  operating  at  a  deviation  angle  of  about  8  degrees.  Lenses  with  focal  lengths 
of  -12.75  and  250  mm  are  used,  in  a  telescope  arrangement,  to  magnify  the  beam  before 
diffraction  by  the  grating.  A  750  mm  lens  focuses  the  diffracted  beam  through  a  50  pm 
slit.  To  obtain  frequency  domain  correlation  curves,  we  fix  the  output  slit  in  Channel  1 
at  Xi  and  scan  the  position  X2  in  Channel  2.  The  photons  that  are  transmitted  through 
the  monochromator  slits  are  coupled  into  multimode  hbers  and  detected  with  time-resolved 
single  photon  counting  modules  (SPCMs,  id  Quantique  id400  and  PerkinElmer  SPCM-AQR- 
16-FC).  A  photograph  of  the  monochromators  is  shown  in  Fig.  18. 

The  primary  experimental  results  of  this  work  are  shown  in  Fig.  19.  For  each  case,  we  set 
the  monochromator  slit  in  Channel  1  at  an  arbitrary  position  xi  that  is  near  the  center  of 
the  generated  32  nm  spectrum  and  leave  the  position  of  this  slit  fixed  thereafter.  The  slit  in 
Channel  2  is  scanned  over  positions  X2,  and  the  coincidence  rate  of  the  two  detectors  (with 
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Figure  18:  Photo  of  frequency  correlation  setup.  Identical  monochromators  measure  the 
frequencies  of  the  entangled  photons.  Red  and  blue  lines  show  the  beam  path  for  Channel  1 
and  Channel  2,  respectively. 

gate  width  T  =  1.25  ns)  is  recorded  as  a  function  of  this  position.  For  each  position,  the 
rate  is  averaged  for  20  s. 

With  the  pump  frequency  dehned  as  ujp,  and  the  position  X2  proportional  to  the  frequency 
UJ2,  we  express  the  coincidence  rate  as  a  function  of  relative  frequency  A  =  uj2  —  —  i^i)- 

The  scale  of  the  frequency  axis  is  calibrated  by  measuring  the  sideband  spacing  of  a  single¬ 
mode  1064  nm  laser  modulated  at  30  GHz,  with  the  zero  position  chosen  (at  the  start  of  the 
experiment)  as  the  location  of  the  correlation  peak  for  unmodulated  photon  pairs. 

The  upper  portion  of  Fig.  19  shows  the  experimental  results  without  modulation  and 
with  modulation  in  a  single  channel.  In  Fig.  19(a)  (upper),  both  modulators  are  turned 
off  by  disconnecting  their  30  GHz  drive  signals.  As  expected  by  energy  conservation,  a 
single  correlation  peak  is  observed.  In  Fig.  19(b)  (upper),  Channel  1  is  phase  modulated 
as  exp[i(5sin(a;mt)]  with  a  modulation  depth  of  5  =  1.5,  and  Channel  2  is  not  modulated. 
The  frequency  correlation  is  now  distributed  over  a  set  of  sidebands,  having  Bessel  function 
amplitudes  J^(5),  whose  total  area  is  equal  to  that  of  upper  (a)  of  Fig.  19. 
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Figure  19:  Frequency  correlation  measurements.  The  horizontal  axes  measure  relative  fre¬ 
quency  changes  as  slit  position  X2  is  scanned.  Upper:  (a)  both  modulators  turned  off  and 
(b)  the  modulator  in  Channel  1  running  at  a  modulation  depth  of  1.5.  Lower:  both  modu¬ 
lators  running  (a)  with  the  same  phase  and  (b)  with  opposite  phase.  The  modulation  depth 
for  both  modulators  is  1.5.  Circles  are  data;  curves  are  theoretical  hts  (see  text). 

In  Fig.  19(a)  (lower),  both  modulators  are  turned  on  at  a  modulation  depth  of  5  =  1.5,  and 
the  cable  length  is  adjusted  so  that  they  have  the  same  phase.  They  now  act  cumulatively 
(constructively  interfere)  to  produce  a  set  of  sidebands  having  a  Bessel  function  distribution 
J^{26).  The  frequency-domain  correlation  function  of  two  distant  modulators  is  therefore 
the  same  as  that  which  would  be  obtained  by  correlating  an  unmodulated  photon  with  a 
photon  modulated  at  twice  the  modulation  depth. 

In  Fig.  19(b)  (lower),  the  modulators  are  run  at  the  same  depth  as  in  the  previous  para¬ 
graph,  but  now  the  relative  cable  length  is  adjusted  so  that  the  modulators  are  run  in  phase 
opposition.  The  modulators  now  destructively  interfere,  and  no  sidebands  are  visible. 

The  striking  difference  between  Figs.  16  and  19  is  that  the  sidebands  are  not  delta  func¬ 
tions  but  have  widths  determined  by  the  resolution  of  the  frequency  correlator  in  Fig.  15 
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(right).  We  have  developed  a  theory  to  predict  the  measured  correlation  functions  with 
hnite- resolution  monochromators  as  used  in  our  experiment.  The  following  describes  this 
theory,  which  is  used  to  produce  the  solid  curves  in  Fig.  19. 

Working  in  the  Heisenberg  picture,  a  nonlinear  crystal  of  length  L  is  pumped  by  a 
monochromatic  laser  at  frequency  ujp.  A  positive-frequency  held  operator  a(a;,  z),  represent¬ 
ing  entangled  photons,  evolves  inside  the  crystal  and  may  be  written  in  terms  of  an  envelope 
b{uj,z)  which  varies  slowly  along  the  propagation  direction:  a{uj,z)  =  b{uj,  z)  exp[ik{uj)z]. 
The  propagation  equations  describing  entangled  photon  generation  are 


db{u,  z) 
dz 

dW{ijj,  z) 
dz 


iK{ijj)b\ujp  —  oj,  z)  exp  [i^k{ijj)z\, 
—ilk*  {uj)b{ujp  —  u,  z)  exp  [—iAk{uj)z] . 


(7.1) 


where  k{uj)  and  Ak{uj)  are  the  coupling  factor  and  wave-vector  mismatch,  respectively.  The 
solution  for  the  output  held  aX  z  =  L,  expressed  in  terms  of  the  vacuum  held  at  the 

input  of  the  crystal,  is 


l^out(^)  A (cU ) fl vac (a.^ )  T  ^(^)^tac(^P  ('^•^) 

where,  to  preserve  the  commutation  relations,  the  functions  A((z)  and  B(u;)  satisfy  \  A[uj)\^  — 
|H(a;)p  =  1  and  A{ijj)B{ijjp  —  uj)  =  B(uj)A{ujp  —  uj). 

The  time-domain  output  held  operator  is  related  to  its  frequency-domain  counterpart 
[Eq.  (7.2)]  by  the  inverse  Fourier  transform,  aout(t)  =  ciout(^)  Gxp(—iLijt)duj,  and  is  nor¬ 
malized  so  that  the  total  rate  of  generated  photons  exiting  the  crystal  is  Rout  =  (olut(^)®out(^))- 
The  generated  photons  are  separated  into  two  channels,  denoted  as  Channel  1  and  Chan¬ 
nel  2,  using  a  50/50  beam  splitter.  The  held  operators  at  the  outputs  of  the  beam  splitter 
are  ai(t)  =  02 (t)  =  :^aout(^)-  The  photons  are  modulated  by  periodic  phase  modula¬ 
tors  whose  time-domain,  Fourier-series  transfer  functions  are  mi{t)  =  qk  exp{—ikumt)  in 
Channel  1  and  m2{t)  =  '^iriexp{—ilumt)  in  Channel  2,  with  Fourier  transforms  mi  (a;)  = 
—  kujm)  and  m2{oj)  =  ~  respectively.  With  the  *  symbol  denot¬ 

ing  convolution,  the  frequency-domain  modulated  helds  are  ai(a;)  =  ai(a;)  *  mi(ci;)  and 
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a2{uj)  =  a2{uj)  *  7712(0;).  Substituting  ai(o;),  02(0;),  7?7i(o;),  and  7772(0;)  into  the  expressions 
for  ai(o;)  and  02(0;)  yields 


di(o;) 


02(0;) 


^  ^  |^y4(o;  kui Yn) OjYajciy^ 


k=—oo 


V2 

+  B{uj  -  kuJm)al^c{ujp  -UJ  +  kUm)]  , 

^  00 

'  /=— 00 

+  B{uj  -  lujrn)al^^{ujp  -  o;  +  luJrn)]  ■ 


(7.3) 


The  modulated  photons  are  frequency  correlated  by  passing  each  through  identical  monochro¬ 
mators  whose  output  slits  may  be  translated  to  select  frequencies  ui  =  (3xi  in  Chan¬ 
nel  1  and  002  =  (3x2  in  Channel  2,  where  the  constant  (3  is  the  linear  dispersion  of  the 
grating  systems.  The  monochromators  (spectral  hlters)  have  held  transmission  functions 
Hi{u}  —  (3xi)  and  H2{uj  —  {3x2)-  The  hltered  held  operators  in  Channels  1  and  2  are 
aif(ct;,a:i)  =  ai(ci;)i7i(ci;  —  (3xi)  and  a2f(ci;,  X2)  =  02(07)7/2(1^  —  (3x2),  respectively.  The  count 
rates  at  the  outputs  of  the  monochromators  are  given  by  i?i(a;i)  =  (ajf(t,  xi)aif(t,  xi))  and 
R2{x2)  =  {a\i{t-,X2)a2i{t,X2)) .  These  rates  are 


2  roo 

Ri{xi)  =  |5(a;  - /ca;m)|^  |77i(o;  - /3a;i)|^(io;, 

k=—oo 

1  noo 

R2{x2)  =  \B{uj -lujm)?\H2{uj  -  Px2)\^du.  (7.4) 

^Z=-oo 

Assuming  a  gate  width  T,  the  coincidence  rate  for  the  two  detectors  is  related  to  the 
second-order  Glauber  correlation  function  G^'^\ti,xi,t2,X2)  = 

(a2f(t2,  a;2)ajf(ti,  a;i)aif(ti,  a;i)a2f(t2, 3;2))-  With  the  assumption  that  the  resolution  of  the 
monochromators  is  high,  or  equivalently  that  the  hlter  widths  are  small  (as  compared  to  the 
modulation  frequency  uJm),  h  can  be  shown  that  the  correlation  function  depends  only  on 
the  difference  of  the  arrival  times  t  =  t2  —  ti,  and  the  coincidence  rate  is 
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Rc{Xi,X2) 


(7.5) 


Equation  (7.5)  may  be  expanded  using  Wick’s  theorem  and  shown  to  be  given  by 


Rc{Xi,X2)  =  Rl{Xi)R2{x2)T 


'^qkTn-kFkir,  xi,  X2) 


k=—oo 


2 

dr, 


(7,6) 


where  A  =  j3{xi  +  X2)  —  cUp,  n  =  [A/uJm  +  |J,  and 


1 

Fk{r,  xi,X2)  =  ^  A{uj  -  kuJm)B{ujp- UJ  +  kujrn) 

Jo 

X  Hi{uj  —  (Jxi)H2{ojp  —  uj  —  j3x2  +  nujm)  exp{iujT)duj.  (7.7) 

The  hrst  term  in  Eq.  (7.6)  is  the  result  of  accidental  coincidences  between  unpaired  pho¬ 
tons  in  a  gate  width  T.  The  second  term  is  the  coincidence  rate  between  paired  photons 
and  captures  the  modulation  effects  described  in  this  section.  To  obtain  Eqs.  (7.4)-(7.7),  we 
have  assumed  that  the  transmission  widths  of  the  monochromators  are  small  as  compared  to 
the  modulation  frequency  and  large  as  compared  to  the  inverse  of  the  temporal  gate  width 
T.  In  our  experiment  these  assumptions  are  satished  by  factors  of  3.5  and  11,  respectively. 

If  we  further  assume  that  A{uj)  and  B[ujp  —  uj)  are  constant  in  the  vicinity  (±150  GHz  in 
Fig.  19)  of  a;  =  {Jxi  and  are  equal  to  Aq  and  Bq,  respectively,  then  Eq.  (7.6)  becomes 


Rc{A)  =  R1R2T  ±  CnH^'^\nujm  -  A), 
where  H^‘^\uj)  =  |i7i(c<;)p  *  |i72(i^)P,  and 


Ri 

R2 

C-n 


1 

Stt 


AqBq 


Hi{uj)\^duj, 

H2{uj)\^duj, 

2 


Qk'^n—k 


k=—oo 


(7,8) 


(7.8a) 

(7.8b) 
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(7.8c) 


The  solid  curves  in  Fig.  19  are  theoretical  hts  to  the  data  using  Eq.  (7.8)  shifted  horizon¬ 
tally  so  as  to  match  center.  The  Fourier  series  coefficients  for  sinusoidal  phase  modulators 
are  Bessel  functions  with  qk  =  Jk{—Si)  and  ri  =  Ji{—S2),  where  (5i  and  S2  are  the  modu¬ 
lation  depths  in  Channels  1  and  2,  respectively  (|(5i|  =  |(52|  =  1-5  in  our  experiment).  We 
model  the  monochromator  response  functions  in  Channels  1  and  2  as  Gaussians  with  FWHM 
bandwidths  F:  Hi{u)  =  ai  exp  [— 2  ln(2)a;^/F^]  and  H2{uj)  =  a!2  exp  [— 2  ln(2)a;^/F^].  (The 
monochromator  in  Channel  1  is  the  mirror  image  of  the  one  in  Channel  2  which  has  a  mea¬ 
sured  FWHM  bandwidth  of  8.5  GHz.)  The  transfer  functions  include  htting  parameters  ai 
and  Q!2  used  in  Fig.  19  to  account  for  transmission  losses  and  the  difference  in  detection 
efficiencies  of  the  photon  counters. 

To  obtain  the  constants  Aq  and  Bq,  for  each  case  in  Fig.  19,  we  measure  the  average 
value  of  i?2  and  use  Eq.  (7.8b)  to  calculate  |i?o|.  We  obtain  |Ho|  from  the  commutator¬ 
preserving  condition  |HoP  —  |5oP  =  1-  For  all  curves,  the  htting  parameters  are  taken  as 
af  =  1.20  X  10“^  and  =  5.59  x  10“'^.  These  values  are  in  good  agreement  with  loss 
measurements  and  estimates  of  the  photon  counter  detection  efficiency,  where  we  note  that 
the  id400  detector  in  Channel  1  has  a  detection  efficiency  an  order  of  magnitude  larger  than 
the  SPCM-AQR-16-FC  detector  in  Channel  2. 

The  experiment  and  theory  summarized  in  this  section  pertain  to  sinusoidal  phase  mod¬ 
ulation.  We  have  also  developed  a  general  theory  for  nonlocal  modulation,  presented  below. 

For  general  modulation,  the  modulated  signal  and  idler  held  operators  are  described  by 
convolution  as  before:  ai(ci;)  =  ai(ci;)  *  mi{uj)  and  02(0;)  =  02(0;)  *  m2{uj);  however,  the 
modulation  functions  mi{t)  and  m2{t)  [with  Fourier  transforms  mi  (a;)  and  m2  (a;)]  are  no 
longer  assumed  to  be  periodic.  The  resulting  Glauber  correlation  function  becomes: 


G^‘^\ti,Xi,t2,X2)  =  C{ti,Xi,t2,X2)  +  Q{tl,Xi,t2,X2), 
where  the  functions  C(ti,xi,t2,X2)  and  Q{ti,xi,t2,X2)  are 


(7.9) 
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^oo  roo  roo  roo  poo  poo 


C{ti,Xi,t2,X2)  = 


IGvr^ 


ml{Qi)mi{Qs)rn2{^2)'iTT‘2{^4) 


'  —CO  J  —CO  J  —CO  J  —CO  J  —CO  J  —CO 

|2  I  D/.  .  O  m2  lj*  / 


Q(ti,a:i,t2,2:2)  = 


X  |i?(a;i  —  fii) I  |i?(a;2~^2)|  (cui  — /3xi)-ffi(a;i  —  fii  +  1^3  — /Jxi) 

XH^  {0O2  -  Px2)H2{uJ2  -  ^2  +  ^4  -  PX2) 

X  exp{z[(f22  —  ^24)^2  +  (^1  —  fi3)ti]}duiduj2dQidQ2dfl3dQ4, 

1 


pco  poo  poo 


7711(^21)7712(^2)^(1^  —  ^i)B{oi!p  —  a;  +  f2i) 


'—00  J —00  J —00 


1677^ 

XHiiiO  —  (3xi)H2{uJp  —  Ci7  +  f2i  +  1^2  —  PX2) 

X  exp[7(a;  —  f2i  —  fl2)it2  —  ti)]  exp[— 7(12i  +  f22)ti](ici;(il2i(il22 


(7.10) 


In  order  to  observe  nonlocal  modulation  interference,  the  modulating  waveform  must  vary 
slowly  as  compared  to  the  temporal  correlation  width  of  the  entangled  photons  (as  described 
by  the  paper  by  Harris  in  2008).  In  Eqs.  (7.10),  the  slow  modulation  condition  may  be 
applied  by  taking  A{uj)  and  B{uj)  to  be  constants  (representing  a  transform-limited,  inhnite- 
bandwidth  biphoton  whose  temporal  correlation  is  a  delta  function).  We  also  assume  that 
the  frequency  correlation  measurement  has  inhnite  resolution:  Hi[uj)  =  H2{uj)  =  6{lj).  With 
these  assumptions,  C{ti,xi,t2,X2)  becomes  a  constant,  independent  of  xi  and  X2  (or,  equiv¬ 
alently,  A);  it  is  a  classical  background  term  arising  from  coincidences  between  uncorrelated 
photons  from  different  entangled  pairs.  Qiti,xi,t2,X2)  becomes 


Q{ti,Xi,t2,X2)  oc  |(7?7i  *7?72)(A)P 

=  \B{mi{t)m2{t)}\‘^ ,  (7.11) 

where  (ttii  *  7772) (a;)  =  7?7i(ci;)  *  7772(147),  and  B{f(t)}  =  ^  f(t)  exp{iujt)dt  is  the  Fourier 
transform.  If  mi{t)  and  7772(f)  are  phase  modulation  waveforms,  i.e.  7771(f)  =  exp[7$i(f)] 
and  7772(f)  =  exp[7<I)2(f)],  then  the  modulators  interfere  so  as  to  produce  a  frequency  domain 
correlation  proportional  to  the  square  of  the  Fourier  transform  of  exp{7[$i(f)  -|-  <l>2(f)]}-  For 
this  relation  to  hold,  it  is  required  that  *Fi(f)  and  <1)2  (f)  both  vary  slowly  as  compared  to  the 
temporal  width  of  the  two-photon  wavefunction. 
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8  Hiding  Single  Photons  with  Spread  Spectrum  Tech¬ 
nology 

In  this  section  we  describe  a  proof-of-principal  experiment  demonstrating  the  use  of  spread 
spectrum  technology  at  the  single  photon  level.  We  show  how  single  photons  with  a  pre¬ 
scribed  temporal  waveform  may  be  transmitted  through  a  noisy  environment  that  is  created 
by  either  narrow  band  thermal  photons,  or  as  in  this  work,  an  interfering  laser  beam  that 
has  an  average  power  that  is  about  a  thousand  times  larger  than  the  average  power  of  the 
beam  of  single  photons.  We  do  this  by  using  two  synchronously  driven  electro-optic  phase 
modulators.  As  shown  in  Fig. 20 (a),  the  modulator  at  the  transmitter.  Mi  broadens  the 
spectrum  of  the  incident  photon  beam  from  about  1  MHz  to  about  10  GHz  and  thereby 
reduces  the  spectral  power  density  by  a  factor  of  10^.  The  receiver  modulator,  M2  is  run  in 
anti-phase  to  the  transmitter  modulator  so  as  to  demodulate  the  photon  beam  and  reduce 
its  bandwidth  to  the  original  1  MHz,  thereby  allowing  it  to  be  transmitted  through  a  narrow 
bandpass  hlter,  F2. 

Now  suppose  that  one  wishes  to  transmit  the  single  photon  beam  through  an  environment 
of  thermal  photons  that  have  a  linewidth  that  is  comparable  to  that  of  the  original  photon 
beam.  If  this  noise,  or  instead,  a  narrowband  interfering  laser  beam  is  injected  after  the  hrst 
modulator  [Fig. 20(a)],  it  is  spectrally  broadened  by  the  second  modulator  and  only  a  fraction 
of  its  power  is  transmitted  through  the  narrow  band  hlter  at  the  receiver.  As  a  result,  the 
signal-to-noise  is  increased  approximately  in  the  ratio  of  the  modulation  bandwidth  to  the 
bandwidth  of  the  hnal  hlter,  and  in  the  language  of  spread  spectrum  technology  the  system 
experiences  a  “processing  gain” .  In  this  work  we  demonstrate  a  processing  gain  of  a  factor 
of  50. 

In  this  experiment  we  use  narrowband  time-energy  entangled  pairs  of  photons  (biphotons) 
produced  as  described  earlier.  We  set  the  optical  depth  of  the  MOT  at  30.  With  2713  equal 
to  the  spontaneous  decay  rate  out  of  state  |3),  we  set  the  coupling  laser  Rabi  frequency 
He  =  4.OI713,  the  pump  laser  Rabi  frequency  Qp  =  1. 66713  and  the  pump  laser  detuning, 
Aujp  =  48.67713.  Under  these  conditions  the  source  produces  narrowband  biphotons  with  an 
estimated  linewidth  of  3.5  MHz  at  a  rate  of  about  11500  pairs/s. 
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Figure  20:  Schematic  of  paired  photon  generation  and  spread  spectrum  experiment,  (a)  The 
spectrum  of  anti-Stokes  photons  produced  in  a  magneto-optic  trap  is  spread  and  narrowed 
by  phase  modulators  Mi  and  M2  driven  by  a  pseudorandom  bit  sequence  (PRBS)  generator. 
A  weak  narrowband  laser  is  used  to  hide  the  photons.  Fi  and  F2  are  narrowband  hlters.  A 
50/50  beamsplitter  (BS)  is  used  to  verify  the  presence  of  single  photons,  (b)  Pseudorandom 
phase  modulation  in  the  time  domain.  The  labels  1,  0  refer  to  the  outputs  of  the  PRBS 
generator,  (c)  Schematic  of  the  optical  spectrum  after  phase  modulation  by  the  PRBS 
generator. 


The  detection  of  the  Stokes  photon  by  detector  Di  sets  the  time  origin  for  measurement 
of  the  biphoton  wavefunction  at  detector  ^2.  The  anti-Stokes  photon  is  sent  through  a 
20-GHz  electro-optic  phase  modulator  (Eospace  Inc.)  driven  by  one  of  two  differential  data 
outputs  of  a  pseudorandom  bit  sequence  (PRBS)  generator  (Centellax  TG2P1A)  with  a  bit 
rate  tunable  between  0.05-10  Gb/s.  In  each  time  slot  AT,  equal  to  the  inverse  bit  rate,  the 
PRBS  outputs  a  random  voltage  equal  to  14/2  (corresponding  to  logic  1)  or  —14/2  (logic 
0),  where  14  is  the  half-wave  voltage  of  the  phase  modulators.  We  choose  a  bit  pattern 
consisting  of  32767  (2^®  —  1)  random  bits.  At  a  bit  rate  of  10  Gb/s,  the  phase  therefore 
switches  randomly  [Fig.20  (b)]  between  — 7r/2  and  7i/2  radians  every  0.1  ns.  This  process  of 
random  phase  modulation  in  the  time  domain  spreads  the  spectrum  in  the  frequency  domain 
from  about  3.5  MHz  for  the  unmodulated  photon  to  about  10  GHz  for  the  modulated  photon. 
The  envelope  of  the  spread  spectrum  is  a  sine-squared  function  with  the  hrst  null,  as  shown 
in  Fig.20  (c),  at  a  frequency  equal  to  the  bit  rate. 
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The  phase-modulated  anti-Stokes  photon  is  transmitted  over  a  channel  consisting  of  a 
long  single-mode  optical  hber  [Fig. 20(a)].  At  the  receiver  end,  the  photon  is  sent  through  a 
second  phase  modulator,  identical  to  the  hrst  and  driven  by  the  complementary  output  of 
the  PRBS  generator.  The  two  phase  modulators  are  synchronized  by  adjusting  the  lengths 
of  the  electric  cables  connecting  them  to  the  PRBS  generator  such  that  the  relative  time 
delay,  Atup  of  the  radio  frequency  signal  sent  to  the  two  modulators  equals  the  travel 
time,  Atp,  of  the  anti-Stokes  photon  between  the  modulators.  Fine  adjustment  of  At^ip  is 
achieved  by  means  of  a  radio  frequency  waveguide  of  variable  length  connected  to  one  of 
the  two  rf  cables.  When  the  phase  modulators  are  run  in  synchronism  such  that  the  phase 
imposed  on  the  photon  by  the  hrst  modulator  is  exactly  undone  by  the  second  modulator, 
the  spectral  width  of  the  photon  is  restored  to  its  original  value  of  3.5  MHz.  When  the 
phase  modulators  are  not  run  in  exact  phase  opposition,  the  uncompensated  residual  phase 
results  in  a  spectrum  that  is  wider  than  3.5  MHz.  After  passing  through  the  second  phase 
modulator  the  anti-Stokes  photon  is  sent  through  a  65-MHz  hber  based  Fabry-Perot  hlter 
(Micron  Optics)  with  a  free  spectral  range  of  13.6  GHz.  If  the  spectral  width  of  the  photon 
after  the  second  phase  modulator  is  less  than  the  hlter  bandwidth,  the  photon  passes  through 
with  little  loss  of  probability  amplitude.  Conversely,  if  the  spectral  width  of  the  photon  is 
larger,  then  it  is  attenuated  by  the  hlter.  The  spectrum  of  narrowband  noise  entering  the 
channel  after  the  hrst  modulator  is  spread  by  the  second  modulator  and  is  attenuated  by  the 
hlter.  To  hide  the  photon  we  use  an  attenuated  beam  from  an  external-cavity  diode  laser 
with  a  linewidth  of  about  300  kHz  set  at  the  same  frequency  as  the  anti-Stokes  photon.  A 
50/50  beamsplitter  (BS)  and  a  third  detector  is  placed  after  the  hlter  to  verify  the  single 
photon  nature  of  the  experiment. 

With  both  phase  modulators  and  the  (noise  simulating)  laser  turned  oh,  the  shape  of  the 
anti-Stokes  photon  wavepacket  is  shown  in  Fig.  21(a).  This  shape  is  measured  by  record¬ 
ing  coincidence  counts  between  detectors  Di  and  D2,  with  the  beamsplitter  removed,  as  a 
function  of  the  diherence  in  arrival  times  of  the  Stokes,  tg,  and  anti-Stokes,  tas,  photons 
at  the  two  detectors.  When  the  hrst  phase  modulator  is  switched  on,  the  spectrum  of  this 
photon  is  spread  to  about  10  GHz.  If  the  second  phase  modulator  is  oh,  the  anti-Stokes 
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Figure  21:  Shape  of  the  anti-Stokes  photon  wavepacket  measured  by  recording  coincidence 
counts  between  detectors  Di  and  D2  with  the  beamsplitter  removed,  (a)  Anti-Stokes  photon 
wavepacket  measured  when  modulators  Mi  and  M2  are  off.  (b) Modulator  Mi  is  on  and 
modulator  M2  is  off.  (c)  Both  modulators  are  on  and  run  synchronously,  (d)  Both  modulators 
are  on  with  >  Atp  (see  text). 

photon  is  severely  attenuated  by  the  narrowband  hlter  F2  and  very  little  photon  probability 
amplitude  leaks  through  to  detector  D2-  This  is  shown  in  Fig. 21  (b).  When  the  second 
phase  modulator  is  switched  on  and  driven  so  as  to  cancel  the  phase  imposed  by  the  hrst 
modulator,  the  photon  spectrum  is  narrowed  and  passes  through  hlter  F2  with  little  loss  of 
probability  amplitude.  This  is  shown  in  Fig. 21(c).  The  small  residual  loss  in  photon  prob¬ 
ability  amplitude  probably  arises  from  a  slight  mismatch  between  the  two  modulators  and 
their  complementary  driving  signals.  In  Fig. 21(d),  the  difference  in  lengths  of  the  electric 
cables  connecting  the  PRBS  to  the  two  modulators  is  chosen  such  that  At^j^  is  greater  than 
Atp  by  about  40ps.  As  a  result,  the  phase  imposed  on  the  photon  by  the  hrst  modulator  is 
not  exactly  negated  by  the  second  modulator.  Now,  after  transmission  through  hlter  F2  we 
observe  a  reduced,  but  not  zero,  coincidence  count  rate.  A  similar  result  is  obtained  when 
AtijF  <  Atp. 

We  next  deliberately  inject  a  weak  diode  laser  beam  to  hide  the  the  anti-Stokes  photon. 
With  additional  laser  photons  at  a  rate  of  40000  and  anti-Stokes  signal  photons  at  a  rate 
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(nanoseconds) 


Figure  22:  Retrieval  of  the  anti-Stokes  photon  from  injected  noise,  (a)  Anti-Stokes  photon 
wavepacket  as  hidden  by  an  additional  laser  beam  that  is  injected  into  the  channel.  Here, 
both  modulators  are  off.  (b)  Switching  on  both  modulators  and  running  them  synchronously 
recovers  the  anti-Stokes  photon 

of  30  [Fig.  22  (a)],  anti-Stokes  photons  can  no  longer  be  seen.  If  the  two  modulators  are 
now  turned  on  and  run  in  synchronism,  the  spectrum  of  the  anti-Stokes  photon  is  spread 
by  the  hrst  modulator  and  compressed  by  the  second  modulator  so  that  the  photon  passes 
through  the  hlter  and  arrives  at  the  detector  with  very  little  attenuation.  The  spectrum  of 
the  laser  photons,  on  the  other  hand,  is  spread  by  the  second  phase  modulator  and  attenuated 
by  the  hlter  and  very  little  of  the  laser  power  reaches  the  detector.  When  this  is  the  case,  the 
shape  of  the  anti-Stokes  photon  wavepacket  can  be  recovered  from  the  noise  [Fig.22(b)].  The 
decrease  in  the  level  of  the  hat  uncorrelated  background  by  a  factor  of  about  50  corresponds 
to  the  spread  spectrum  processing  gain. 

The  quality  of  a  single  photon  source  can  be  quantihed  by  means  of  the  conditional 
Glauber  correlation  function,  5'con(i(0)  dehned  in  Eq  5.2.  We  turn  oh  the  the  weak  narrowband 
laser  beam  and  measure  fi'cond(d)  ^  function  of  the  rate  of  detected  Stokes  photons  with 
both  modulators  turned  oh  (blue  line  in  Fig.  23)  and  with  both  modulators  turned  on  and 
running  synchronously  (red  line  in  Fig.  23).  At  a  Stokes  rate  of  approximately  1200s“^  we 
obtain  gfJndi^)  =  0-22  ±  0.06  for  the  unmodulated  photon  and  fi'coL(O)  =  0.14  ±  0.04  for  the 
phase  modulated-demodulated  photon.  These  values  of  than  the  classical 

limit  of  1  (dashed  line  in  Fig.  23)  and  reasonably  less  than  the  limiting  value  of  0.5;  i.e.,  the 
conditional  Glauber  correlation  function  for  a  two-photon  Fock  state. 

We  have  shown  how  narrowband  single  photons  may  be  phase  modulated  so  as  to  increase 
the  width  of  their  spectrum  by  several  orders  of  magnitude  while  at  the  same  time  retaining 
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Figure  23:  Conditional  threefold  correlation  function  as  a  function  of  Stokes  rate  with  both 
modulators  off  (blue)  and  with  both  modulators  on  (red).  The  dashed  line  is  the  classical 
limit.  The  errors  are  calculated  by  assuming  a  Poisson  distribution  for  the  coincidence  rates. 

the  information  that  characterizes  their  waveform.  By  compensation  with  an  anti-phased 
modulator,  the  biphoton  waveform  may  be  reconstructed  at  a  distant  location.  We  have 
also  shown  how  spread  spectrum  technology  may  be  used  to  hide  a  single  photon  in  the 
presence  of  laser  photons  of  the  same  frequency  and  similar  linewidth.  Applications  may 
include  an  additional  level  of  classical  security  for  quantum  key  distribution,  and  someday, 
multiplexing  of  single  photon  communication  channels.  This  work  may  also  be  the  hrst  use 
of  spread  spectrum  techniques  combined  with  single  photon  detection. 
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B 


Construction  and  Operation  of  the  2-D  MOT 


Figure  24:  Magnetic  field  coil  and  field  geometry  for  a  two-dimensional  magneto-optic  trap 
(2D  MOT),  (a)  The  coil  used  to  produce  the  linearly  varying  magnetic  held,  (b)  Magnetic 
held  in  the  x  —  y  plane,  (c)  Magnetic  held  as  a  function  of  x  coordinate  for  coil  current  of 
62A.  (d)  Gradient  of  the  magnetic  held  as  a  function  of  x  coordinate.  The  held  gradient  at 
the  center  of  the  trap  is  about  10  gauss/cm. 

Trapping  of  neutral  atoms  in  a  MOT  is  accomplished  by  applying  a  magnetic  held  which 
varies  linearly  with  the  spatial  coordinate  around  the  region  where  the  atoms  are  to  be 
trapped.  To  trap  the  atoms  in  a  cigar  shaped  MOT,  we  provide  spatial  conhnement  along 
the  two  transverse  directions  only,  with  no  magnetic  held  gradient  along  the  longitudinal 
axis.  Instead  of  the  anti-Helmholtz  coils  of  a  3D  MOT,  we  use  an  elongated  coil  as  shown  in 
Fig.24(a).  This  coil  provides  a  magnetic  held  gradient  of  about  10  gauss/cm  (at  a  coil  current 
of  62  A)  in  the  transverse  direction  and  an  almost  negligible  gradient  along  the  longitudinal 
2;  direction.  The  coil  is  shaped  like  a  racetrack  and  is  made  with  a  square  copper  wire  of  5 
mm  size  with  a  hollow  core.  Water  is  circulated  through  this  core  to  cool  the  coil  and  allow 
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operation  up  to  about  lOOA.  The  shape  of  the  coil  is  chosen  to  reduce  its  inductance  and  to 
therefore  facilitate  fast  switch  on  and  switch  off  of  the  coil  current.  Fig.  24(b)  shows  a  held 
plot  of  the  magnetic  held  in  the  transverse  plane.  At  the  center  of  the  trap,  the  magnetic  held 
is  zero  and  away  from  this  center  it  increases  linearly  as  a  function  of  the  radial  coordinate. 
Fig.  24  (c)  and  (d)  show  the  magnitude  of  the  magnetic  held  and  its  gradient  as  a  function 
of  the  distance  from  the  center  of  the  trap. 


Figure  25:  Optical  layout  of  2D  MOT  and  biphoton  generation  experiment. 


C  Optical  Layout 

Fig.  25  shows  the  layout  of  the  MOT  beams  and  the  FIT  lasers.  The  trapping  laser  is 
divided  into  four  beams  two  of  which  counterpropagate  along  an  axis  perpendicular  to  the 
axis  of  the  MOT.  The  remaining  two  beams  are  sent  long  an  axis  set  at  45  degrees  to  the 
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MOT  axis.  The  power  distribution  between  the  beams  is  set  to  optimize  the  shape  of  the 
MOT  and  to  ensure  that  the  atom  cloud  expands  isotropically  when  the  magnetic  held  is 
switched  off.  Of  importance,  there  are  no  trapping  beams  and  no  magnetic  held  gradient 
along  the  axis  of  the  MOT.  The  repump  laser  beam  is  sent  at  right  angles  to  the  MOT 
axis.  The  trapping  and  MOT  beams  are  2  cm  in  diameter.  The  pump  laser  beam  has  a 
1/e^  diameter  of  1.46  mm  and  is  circularly  polarized  (cr— ).  The  coupling  laser  beam  has  a 
diameter  of  1.63  mm  and  has  the  opposite  polarization  (cr+).  The  pump  and  coupling  laser 
axis  forms  an  angle  of  2°  with  the  MOT  axis.  Stokes  and  anti-Stokes  photons  are  collected 
along  the  axis  of  the  MOT.  The  Stokes  and  anti-Stokes  hber-to-hber  coupling  efficiency  is 
70%  and  their  1/e^  diameter  is  220  pm.  Quarter- wave  plates  and  polarization  beam-splitters 
are  used  to  select  particular  polarizations  of  the  Stokes  and  anti-Stokes  photons  and  also 
provide  an  additional  level  of  hltering  of  the  generated  photons  from  the  strong  pump  and 
coupling  beams.  To  measure  the  optical  depth  and  EIT,  a  weak  probe  beam  is  sent  through 
the  Stokes  hber  and  its  transmission  though  the  MOT  as  a  function  of  the  detuning  from 
line  center  is  recorded  using  a  photomultiplier  tube  (PMT). 

D  Lasers 

The  wavelengths  of  the  lasers  used  in  MOT  preparation  depends  upon  the  choice  of  the 
atoms  used  in  the  trap.  Our  2D  MOT,  as  originally  designed,  trapped  and  some  of 
the  early  experiments  were  performed  on  this  system.  In  order  to  implement  an  optical 
pumping  scheme,  we  later  switched  to  a  MOT  and  some  later  experiments  described 
in  this  report  have  been  done  with  atoms.  Rubidium  was  chosen  because  of  the  large 
hyperhne  splitting  of  the  5Pi/2  and  5P3/2  levels.  In  alkali  atoms  of  lower  atomic  weights,  the 
hyperhne  splitting  is  of  smaller  magnitude.  Another  important  factor  that  led  to  the  choice 
of  Rnbidinm  was  the  ready  availability  of  commercial  lasers  at  the  Rubidium  Di  and  D2 
transitions. 

Our  optical  system  consists  of  several  lasers  with  ontput  frequencies  locked  to  various 
Rubidium  transitions.  For  laser  cooling  and  trapping  and  for  biphoton  generation  it  is 
of  utmost  importance  that  the  frequencies  of  these  lasers  be  controlled  to  a  high  degree  of 
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Figure  26:  Locking  mechanism  of  the  MOT  and  FIT  lasers.  A  reference  laser  (right)  is  hrst 
locked  to  a  Doppler  free  peak  of  using  a  saturation  spectroscopy  set-up.  Light  from  the 
reference  laser  is  mixed  with  a  portion  of  the  light  coming  from  the  science  laser  (left)  on  a 
fast  photodetector  to  generate  a  beat  note.  The  beat  note  is  stabilized  using  a  phase-locked 
loop  to  lock  the  science  laser  relative  to  the  reference  laser. 

precision.  In  the  experiment  this  is  done  by  locking  hrst,  a  reference  laser  to  one  of  the  peaks 
of  a  Doppler  free  saturation  spectroscopy  spectrum.  The  other  lasers  are  locked  relative  to 
this  reference  laser  by  mixing  a  part  of  the  output  of  these  lasers  with  the  reference  laser 
beam  and  then  using  phase-locked  loops  to  stabilize  the  resulting  beat  notes. 

A  schematic  of  the  locking  scheme  is  shown  in  Fig.  26.  A  small  portion  of  the  light 
coming  from  an  external  cavity  diode  laser,  referred  to  henceforth  as  the  reference  laser,  is 
sent  through  an  electro-optic  phase  modulator  (Thorlabs  EO-PM-NR-Cl,  half-wave  voltage 
lA  ~  175V)  driven  with  a  sinusoidal  voltage  with  a  frequency  of  100  MHz  and  a  peak-to-peak 
amplitude  of  about  50V.  As  a  result  of  this  phase  modulation,  the  instantaneous  frequency 
of  the  laser  light  varies  sinusoidally  as  a  function  of  time.  This  beam  is  then  sent  through  a 
room  temperature  Rubidium  cell  to  a  photodiode  (PDl).  To  see  the  Doppler-free  saturation 
spectroscopy  peaks,  a  portion  of  this  phase  modulated  beam  is  reflected  back  into  the  cell 
as  shown  in  the  hgure.  For  lasers  locked  to  the  Rubidium  D1  line  (795nm),  the  reference 
laser  is  locked  to  the  crossover  peak  of  the  F  =  3  ^  F'  =  2  and  F  =  3  ^  F'  =  3  transitions 
of  ^^Rb.  For  the  D2  line  (780nm),  the  reference  laser  is  locked  to  the  crossover  peak  of  the 
F  =  3  ^  F'  =  2  and  F  =  3  ^  F'  =  4  transitions  of  ^^Rb.  The  lasers  are  locked  to  the 
maxima  of  the  crossover  peaks  using  the  technique  of  Pound,  Drever  and  Hall. 

After  the  reference  laser  has  been  locked,  a  portion  of  the  light  from  the  laser  is  mixed 
with  a  part  of  the  beam  of  the  science  laser  as  shown  in  Fig.  26.  The  photodiode  PD2 
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Figure  27;  Lasers  used  to  prepare  MOT.  The  trapping  laser  cools  and  traps  atoms.  The 
repump  laser  pumps  atomic  population  to  the  upper  ground  state  for  laser  cooling,  (a)  Rb 
energy  level  diagram  for  D2  line  and  laser  positions.  The  trapping  laser  is  red-detuned  by  20 
MHz  from  the  cycling  |2)  — >  |3)  transition.  The  repump  laser  keeps  the  F  =  1  state  empty, 
(b)®^  Rb  energy  level  diagram  for  D2  line  and  laser  positions  for  cooling  and  trapping. 

outputs  a  beat  note  at  a  frequency  equal  to  the  difference  between  the  frequencies  of  the 
reference  and  science  lasers.  This  beat  note  is  then  frequency  stabilized  to  that  of  a  local 
oscillator  using  a  phase-locked  loop.  The  frequency  of  the  local  oscillator  can  be  varied  to 
tune  the  science  laser  to  the  desired  Rubidium  transition. 

The  primary  laser  used  for  laser  cooling  is  a  Titanium  Sapphire  laser  (Coherent  MBR 
110)  pumped  by  a  diode  pumped  ND-YAG  laser  (Coherent  Verdi  VIO).  This  laser  has  a 
peak  output  of  about  650  mW  and  is  red  detuned  by  about  20  MHz  from  the  |5S'i/2R  = 
2)  — >  I5P3/2,  F'  =  3)  transition  of  to  trap  atoms  of  this  isotope.  Due  to  the  proximity 
of  the  other  levels  of  the  5P3/2  manifold,  there  is  a  hnite  probability  of  exciting  the  atoms 
to  the  P'  =  1  and  P'  =  2  levels.  Atoms  thus  excited  may,  through  spontaneous  emission, 
decay  down  to  the  |5S'i/2,  P  =  1)  level  and  are  then  no  longer  cooled  by  the  trapping  laser. 
To  return  these  atoms  back  to  the|5S'i/2,P  =  2)  state,  a  second  laser,  called  the  repump 
laser  tuned  on  resonance  with  the  |5S'i/2P  =  1)  — >  I5P3/2,  F'  =  2)  transition  is  used.  Fig.  27 
shows  the  positions  of  the  trapping  laser  and  repump  lasers  used  to  trap  (a)  and 
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Figure  28:  Locking  position  of  the  pump  laser  used  in  biphoton  generation  for  (a)  and 
(b). 

(b)  atoms.  A  780nm  external  cavity  diode  laser  (ECLD)  (Sacher  Lasertechnik  Littrow  TEC 
series)  with  a  peak  power  of  ~  200  mW  is  used  as  the  repump  laser.  These  lasers  are  locked 
using  the  technique  described  earlier. 

In  Fig.  28  and  Fig.  29,  we  show  the  locking  positions  of  the  lasers  used  for  EIT  mea¬ 
surements  and  for  biphoton  generation.  The  pump  laser  (Sacher  Lasertechnik  Littman  TEC 
series)  operates  at  780nm  and  is  blue-detuned  146  MHz  from  the  5S'i/2F"  =  2  — >  5P3/2F'  =  3 
transition  of  or  the  5S1/2F  =  1  — »■  5P3/2F'  =  2  transition  of  as  shown  in  Fig. 
28.  The  coupling  laser  (New  Focus  Vortex  diode  laser)  operating  at  795nm  is  locked  on 
the  D1  line  as  shown  in  Fig.  29.  The  probe  laser  beam  is  generated  by  sending  a  portion 
of  the  coupling  laser  light  through  a  double-pass  acousto-optic  modulator  (AOM)  operated 
at  a  frequency  of  ~  1.57  GHz  (for  ^^Rb)  or  3.47  GHz  for  Kb.  The  probe  frequency 
can  be  swept  by  40  MHz  on  either  side  of  the  transition  to  perform  optical  depth  and  EIT 
measurements. 
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(a)  (b) 


Figure  29:  Locking  position  of  the  coupling  and  probe  lasers  used  in  FIT  measurements  and 
biphoton  generation  for  (a)  and  (b). 

E  Vacuum  System 

The  MOT  is  prepared  in  a  ultra-high  vacuum  system  consisting  of  an  octagonal  ceramic  and 
glass  cell  [Fig.  30  (a)]  mounted  on  a  six- way  cross  connected  to  a  pumping  station.  The 
pumping  system  consists  of  a  roughing  pump,  a  turbo  pump  and  an  ion  pump.  The  vacuum 
assembly  is  shown  in  Fig.  30(b).  The  turbo  pump  (Varian  Turbo- V70)  is  used  during  initial 
pump  down  and  when  baking  the  system.  It  is  then  switched  off  and  the  ion  pump  (Varian 
Vaclon  Plus  40  StarCell)  is  used  to  obtain  a  vacuum  of  about  10“®  torr.  Small  leaks  in 
the  interface  between  the  windows  and  the  alumina  cell  limit  the  vacuum  to  this  value.  In 
the  absence  of  such  leaks,  vacua  of  the  order  of  10“^°  should  be  achievable.  To  one  port  of 
the  six-way  cross  an  electrical  feed-through  is  connected  and  four  Rubidium  dispensers  are 
spot  welded  between  the  terminals.  When  a  current  is  sent  through  these  terminals,  the 
dispensers  are  heated  to  a  high  temperature  and  release  an  isotopic  mixture  of  and 
^’^Rb  atoms  into  the  vacuum  system.  Depending  on  the  positions  of  the  laser  locks,  either  of 
these  atomic  species  can  be  trapped  in  the  MOT. 


49 


Figure  30:  (a)  Octagonal  vacuum  cell.  The  cell  is  machined  from  alumina  and  anti-reflection 
coated  windows  are  bonded  onto  its  surfaces,  (b)  The  vacuum  assembly  showing  the  oc¬ 
tagonal  cell  and  vacuum  pumps.  Rubidium  dispensers  are  spot-welded  to  an  electrical  feed¬ 
through  connected  to  one  flange  of  the  six-way  cross. 

F  Experimental  Cycle 

The  lasers  locked  to  the  various  Rubidium  lines  are  switched  on  and  off  in  a  hxed  sequence 
during  the  biphoton  generation  experiment.  The  switching  of  the  laser  beams  is  done  using 
AOMs  (ISOMET  1206C)  operating  at  110  MHz  and  with  a  transmission  of  ~  70-80%  in 
the  hrst  order  diffracted  beam.  The  AOMs  are  driven  by  home-built  RF  drivers  and  the 
relative  timing  of  the  lasers  is  controlled  by  a  series  of  timing  boxes  (Stanford  Research 
Systems  DG535  Digital  Pulse  Generator)  which  have  a  resolution  of  1  ps.  The  duration  of 
the  experiment  is  divided  into  time  slots  of  length  equal  to  the  inverse  of  what  we  call  the 
repetition  rate.  When  the  repetition  rate  is  200  Hz,  for  example,  the  experimental  cycle 
i.e.  the  length  of  this  time  slot  is  5  ms.  During  this  5  ms  cycle,  the  trapping  and  repump 
lasers  are  kept  on  for  about  4.5  ms  to  prepare  the  MOT.  During  the  remaining  500  ps,  the 
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Figure  31:  Timing  sequence  of  laser  and  detector  operation.  The  cooling  lasers  are  switched 
off  during  the  experimental  window.  All  times  are  in  microseconds. 

trapping  and  repump  lasers  are  switched  off  and  the  pump  and  coupling  lasers  are  switched 
on  for  the  biphoton  generation  experiment  or  the  coupling  and  probe  lasers  are  turned  on 
for  EIT/OD  measurements.  For  all  the  experiments  described  in  this  report,  the  MOT 
magnetic  held  continued  to  remain  on  through  the  entire  experimental  cycle.  A  schematic 
of  the  experimental  cycle  is  shown  in  Fig. 31.  The  cycle  begins  with  the  arrival  of  a  trigger 
signal  at  t  =  0.  The  trapping  and  repump  lasers  are  switched  off  150  /rs  after  the  arrival 
of  the  trigger  pulse.  A  depump  beam  is  then  turned  on  for  150  /is  to  transfer  the  atomic 
population  from  the  upper  5S'i/2  state  to  the  lower  551/2  state.  The  coupling  laser  and  pump 
lasers  are  then  switched  on  0.5  /is  and  10  fis  respectively  after  the  depump  beam  has  been 
turned  off.  The  cloud  of  Rubidium  atoms  now  starts  emitting  paired  Stokes  and  anti-Stokes 
photons.  The  single  photon  counting  modules  (SPCM)  and  the  time-to-digital  converter 
(TDC)  are  then  activated  to  record  the  coincidences  between  the  photons  as  shown  in  the 
hgure.  The  photon  counting  equipment  and  the  pump  and  coupling  lasers  are  turned  off 
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Figure  32:  The  two-dimensional  MOT  atom  cloud  at  left  and  the  atom  cloud  imaged  through 
the  vacuum  cell  on  the  right. 

and  the  cooling  lasers  are  turned  back  on  650  /rs  after  the  arrival  of  the  trigger  pulse  and 
continue  to  remain  on  till  the  end  of  the  experimental  cycle.  In  this  manner,  the  MOT  is 
operated  with  a  10  percent  duty  cycle  with  the  MOT-preparation  time  equal  to  4.5  ms  and 
the  experimental  window  for  photon  counting  equal  to  about  500  ps.  To  measure  the  optical 
depth  and  FIT  curves,  probe  and  coupling  lasers  are  used  during  the  experimental  window. 
Their  timing  is  shown  in  the  hgure.  The  other  lasers  and  the  photon  counting  equipment 
remain  off  during  these  measurements. 


G  MOT  Characteristics 

A  picture  of  the  2D  MOT  is  shown  in  Fig.  32.  The  cloud  is  about  1.7  cm  in  length  and 
has  an  aspect  ratio  of  25.  The  optical  depth  of  the  MOT  increases  as  the  current  to  the 
Rubidium  getters  is  increased.  The  maximum  optical  depth  observed  in  our  experiments 
is  62.  From  the  optical  depth  measurement,  we  estimate  that  there  are  approximately  10® 
atoms  trapped  in  the  MOT. 

An  optical  depth  measurement  on  the  |5S'i/2,R  =  2)  — >  I5P3/2,  R'  =  3)  transition  of 
is  shown  in  Fig.  33  (a).  The  blue  trace  is  the  transmission  of  the  probe  beam  as  a 
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Figure  33:  Optical  depth  and  FIT  measurements  of  the  2D  MOT.  (a)  Transmission  of  the 
probe  beam  as  a  function  of  detuning  from  line  center  with  (red)  and  without  (blue)  the 
coupling  laser,  (b)  Optical  depth  as  a  function  of  the  dispenser  current. 

function  of  detuning  from  line  center  (0  MHz)  without  the  coupling  laser.  The  optical  depth 
as  estimated  by  a  nonlinear  curve  £t  to  this  data  is  45.  When  coupling  laser  is  switched  on, 
a  narrow  transparency  window  opens  up  at  line  center  as  shown  by  the  red  trace.  The  width 
of  the  transparency  window  is  about  1.2  MHz  which  is  smaller  than  the  natural  linewidth 
of  Rubidium  (5.9  MHz).  We  thus  have  subnatural  FIT.  In  Fig.  33  (b)  we  show  the  optical 
depth  measured  on  this  transition  as  a  function  of  the  current  sent  through  the  Rubidium 
dispenser.  To  reduce  the  risk  of  exhausting  the  Rubidium  dispensers,  we  always  operate  the 
dispensers  below  5A. 
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Chapter  4 

Theory  of  EIT  based  Paired 
Photon  Generation 


Correlated  and  entangled  photon  pairs  are  widely  nsed  in  qnantum  commnnication, 
qnantnm  cryptography  [31]  and  qnantnm  imaging  [32,  33].  These  photon  pairs  are 
usnally  prodnced  via  spontaneous  parametric  down  conversion  in  nonlinear  crystals. 
A  few  years  ago  a  new  approach  to  the  generation  of  paired  photons  has  been  exper¬ 
imentally  demonstrated  by  two  groups,  both  use  electromagnetically  induced  trans¬ 
parency  (EIT)  to  generated  paired  photons  in  an  otherwise  opaque  atomic  medium. 
Working  with  hot  atoms,  Lukin  and  colleagues  have  demonstrated  correlation  be¬ 
tween  generated  pulses  of  light,  as  well  as  storage  and  delayed  extraction  [34].  Work¬ 
ing  with  a  MOT,  Kimble  and  colleagues  have  shown  the  generation  of  nonclassical 
photon  pairs  with  a  programmable  delay  [35] .  Recently,  the  Harris  research  group  at 
Stanford  has  demonstrated  generation  and  rudimentary  waveform  control  of  narrow 
band  biphotons  [36].  More  recently,  Kolchin  and  colleagues  have  shown  paired  pho¬ 
ton  generation  with  a  single  pump  beam  in  a  right  angle  geometry  and  paired  photon 
generation  in  the  ensemble  of  cold  two-level  atoms  [37] .  A  long  coherence  length  and 
a  controllable  bandwidth  of  the  generated  paired  photons  are  the  advantages  of  the 
new  approach,  which  might  be  useful  for  such  applications  as  long  distance  quantum 
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communication  [38]  as  well  as  biplioton  waveform  control  and  shaping  ^ 


4.1  Theory  of  paired  photon  generation 

This  chapter  describes  the  theory  of  correlated  paired  photon  generation  in  a  collec¬ 
tion  of  double  A-type  atoms  [39].  Using  the  Heisenberg-Langevin  method  we  evaluate 
and  analyze  spectral  characteristics  of  the  generated  Stokes  and  anti-Stokes  photons 
and  their  time-correlation  properties.  Of  importance,  we  predict  the  regime  when 
the  correlation  time  of  generated  photons  is  determined  by  the  group  delay  caused 
by  EIT.  This  theoretical  prediction  as  well  as  some  theoretical  aspects  of  this  work 
has  been  outlined  in  Ref.  [36].  The  complete  theoretical  treatment  of  paired  photon 
generation  with  and  without  the  approximation,  that  the  atomic  population  remains 
in  the  ground  state,  has  been  given  in  Ref.  [39]. 


Figure  4.1:  Energy  level  diagram  and  schematic  for  spontaneous  backward- wave 
paired  photon  generation  in  an  atomic  cloud  formed  by  double-A  type  atoms.  In 
the  presence  of  the  pump  and  coupling  lasers  phase-matched,  counter-propagating 
Stokes  and  anti-Stokes  photons  are  generated  into  opposite  directions. 

The  schematic  of  the  process  considered  here  is  shown  in  Fig.  5.7.  In  the  presence 
two  cw  beams  termed  as  the  pump  and  coupling  lasers  with  frequencies  Up  and  cUc, 

^In  a  right-angle  geometry,  where  Stokes  and  anti-Stokes  photons  are  generated  and  collected 
at  right  angles  from  the  direction  of  the  pump-coupling  axis,  applying  the  absorption  mask  on  the 
pump  beam  allows  us  to  create  a  biphoton  with  a  precribed  waveform. 
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paired  spontaneous  photons  termed  as  Stokes  and  anti-Stokes  are  generated  in  the 
atomic  cloud  and  propagate  in  opposite  directions  along  the  z  axis.  In  order  to  keep 
the  parametric  gain  small  we  choose  the  pump  beam  to  be  weak  and  detuned  from  the 
resonance  transition  |1)  — >■  |4).  The  intense  coupling  beam  is  tuned  to  resonance  with 
the  |2)  — >■  |3)  transition  to  enhance  the  atom-held  interaction  and  provide  EIT  for  the 
generated  anti-Stokes  photon.  Under  such  conditions  we  expect  a  small  fraction  of  the 
atomic  population  to  be  in  the  excited  states  while  most  of  it  remains  in  the  ground 
state  |1).  The  frequencies  of  the  generated  photons  obey  the  energy  conservation 
+  ^as  =  ^p  +  ^c-  In  fhe  presence  of  EIT  the  anti-Stokes  photon  escapes  out  of  the 
atomic  cloud  with  very  slow  group  velocity. 

We  note  the  connection  to  earlier  work:  Two-photon  entanglement  in  type-II 
SPEC  has  been  analyzed  [40].  The  possibility  of  quantum  correlated  and  squeezed 
helds  in  backward  wave  EIT  system  has  been  predicted  [41]  and  large  parametric 
gain  and  oscillations  have  been  observed  [42].  Control  of  single  photons  has  been 
discussed  [43].  We  also  note  the  early  studies  on  double-A  atomic  systems  [44]. 

In  this  chapter  we  review  and  discuss  in  detail  the  system  dynamics  putting  em¬ 
phasis  on  the  influence  of  the  EIT  window  on  paired  photon  generation  bandwidth. 
We  show  that  at  low  parametric  gain  the  atomic  system  can  operate  in  two  different 
regimes.  In  the  hrst  regime,  where  the  group  delay  is  small,  the  intensity  correlation 
function  shows  Rabi  oscillations.  In  the  second  regime,  where  the  optical  depth  of  the 
atomic  sample  and  the  group  delay  are  large,  phase-matching  becomes  the  dominant 
process  that  controls  the  shape  of  the  intensity  correlation  function.  We  examine 
the  conditions  required  for  the  system  to  operate  in  the  oscillatory  and  group  delay 
regimes  and  discuss  corresponding  Stokes  and  anti-Stokes  spectral  generation  rates. 
In  particularly,  we  predict  that  the  group  delay  regime  requires  the  optical  depth 
much  higher  that  10.  This  point  has  been  missed  in  our  earlier  publication  [36]. 

We  also  extend  our  theoretical  treatment  of  paired  photon  generation  to  go  beyond 
the  ground  state  approximation.  This  allows  us  to  properly  include  and  analyze  the 
effect  of  Langevin  noise  fluctuation  on  the  atomic  system  and  solve  the  problem  of  its 
return  to  the  ground  state  after  the  emission  of  Stokes  and  anti-Stokes  photon  pairs. 
We  also  address  the  important  questions:  1)  ‘What  are  the  effects  of  the  optical 
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thickness  of  the  atomic  sample  on  paired  and  single  photon  generation?”  2)  “Does 
every  Stokes  photon  have  its  paired  anti-Stokes  photon?” 


We  introduce  and  derive  the  biphoton  wavefunction  taking  into  account  Langevin 


noise  terms.  We  also  derive  and  discuss  the  relation  of  the  Stokes  photon  count  rate 
to  the  classical  gain  coefficient  in  the  Stokes  channel. 

4.2  Heisenberg-Langevin  description  of  paired  pho¬ 
ton  generation 

We  consider  a  collection  of  identical  double  A-type  atoms  uniformly  distributed  within 
a  pencil-shaped  volume  with  cross  section  S  and  length  L.  We  assume  that  the  atomic 
sample  is  optically  thin  in  the  transverse  direction,  so  that  there  is  no  radiation 
trapping  effect  in  this  direction.  No  restrictions  are  imposed  on  the  optical  thickness 
of  the  atomic  sample  in  the  z  direction.  We  also  assume  that  the  pump  and  the 
coupling  beams  counter-propagate  undepleted  through  the  atomic  medium.  Under 
these  assumptions  we  consider  propagation  of  a  single  transverse  spatial  mode  of 
radiation  along  the  z  axis.  The  pump  and  coupling  laser  beams  are  treated  as  classical 
quantities  and  their  interaction  with  the  medium  is  described  semi-classically.  In  order 
to  allow  for  the  spontaneous  initiation  of  the  parametric  fluorescence  process,  the 
generated  weak  Stokes  and  anti-Stokes  helds  are  described  by  quantum-mechanical 
operators,  in  slowly  varying  envelope  approximation: 


(4.1) 


where  subscript  j  denotes  either  Stokes  or  anti-Stokes  photon,  Us  =  —  UJ2  +  Aa;i4, 

=  1^3  —  a;i,  U  =  L  X  S'  is  the  interaction  volume. 

Adopting  the  notations  of  Lukin  and  Fleischhauer  [45,  46],  in  the  rotating  wave 
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approximation  we  write  the  interaction  Hamiltonian  in  continuous  form  as 


1/  = 


hN 


dzi  Auji4a44{z,t)  +  g3idas{z,t)a3i{z,t) 


+  yO-32(^,  t)  +  t)  +  942ds{z,  t)a42{z,  t)  +  h.c.  ) ,  (4.2) 


where  ajk{z,t)  are  the  collective  slowly  varying  atomic  operators,  dehned  in  the  Ap¬ 
pendix  4.7.1,  N  is  the  total  number  of  atoms  in  the  atomic  ensemble,  Aa;i4  is  the 
detuning  of  the  pump  laser  from  the  |1)  — >•  |4)  transition,  Dc  = 
are  pump  and  coupling  laser  Rabi  frequencies  with  Ep,  Ec  as  the  complex  amplitudes 
of  the  electric  helds,  Qs  =  and  Qas  =  are  the  coupling  constants  with  pjk 

as  the  dipole  moment  for  the  \j)  — >•  \k)  transition  and  Sj  =  \J~^^  as  the  electric  held 
of  a  single  photon. 

The  propagation  of  the  Stokes  and  anti-Stokes  helds  and  their  interaction  with 
the  atoms  are  described  by  the  set  of  Maxwell  and  Heisenberg-Langevin  equations. 
The  Heisenberg-Langevin  equations  are  responsible  for  the  atomic  evolution: 


d  . 


+  F,, 


jk-, 


(4.3) 


where  7^^  are  the  dephasing  rates,  are  the  spontaneous  emission  rates,  Fjk{z,t) 
are  the  collective  atomic  ^-correlated  Langevin  noise  operators.  The  full  set  of  the 
Heisenberg-Langevin  equations  is  shown  explicitly  in  Sec.  4.7.2 

The  huctuations  of  (5-correlated  collective  Langevin  noise  operators  Fjk{z,t)  is 
given  by 

{Fjk{z,t)Fyk'{z',t'))  =  ^'Djk,j'k'{z,t)6{t  -  t')S{z  -  z'),  (4.4) 

where  Djkj'k'  is  a  Langevin  dihusion  coefficient.  The  derivation  of  Eq.(4.4)  and 
relevant  dihusion  coefficients  is  shown  in  Sec.  4.7.4. 
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The  evolution  of  the  annihilation  and  creation  operators  for  the  slowly  vary¬ 
ing  Stokes  and  anti-Stokes  helds  is  described  by  the  coupled  propagation  equations 


m 

m 


“4) 


ds{z,t)  =  igsNa2i{z,t) 

=  -igasNdsiizR). 


(4.5) 


The  following  analysis  of  the  system  involves  calculation  of  the  expected  values 
of  quantum  held  operators  and  their  combinations.  We  note  that  in  the  Heisenberg 
picture  operators  evolve  and  the  system  is  always  in  its  initial  state,  which  in  our  case 
corresponds  to  no  Stokes  and  anti-Stokes  input  beams  at  the  left  and  right  boundaries. 


4.3  Solution  of  Coupled  Equations 

Due  to  complexity  and  nonlinearity  it  is  not  possible  to  obtain  an  analytic  solution 
for  the  combined  set  of  the  held  propagation  equations  and  Heisenberg-Langevin 
equations.  Nevertheless,  under  the  condition  that  Stokes  and  anti-Stokes  helds  are 
much  weaker  than  the  coupling  and  pump  helds,  and  Stokes  and  anti-Stokes  photon 
densities  are  much  smaller  than  the  atomic  density  N/V  [45],  the  Heisenberg-Langevin 
equations  can  be  linearized  with  Stokes  a*  and  anti-Stokes  das  helds  as  perturbation 
parameters.  The  linearization  procedure  is  described  in  detail  in  the  Appendix  4.7.3. 

In  order  to  solve  the  set  of  the  linearized  Heisenberg-Langevin  Eq.  (4.35)  and 
coupled  propagation  Eq.  (4.5),  we  hrst  Fourier  transform  them.  Then,  extracting  the 
solutions  for  5'24(a;),  d3i(a;)  and  substituting  them  into  Fourier  transformed  Eq.(4.5), 
we  obtain  the  conpled  eqnations  for  ds{z,uj)  and  dlg{z,  —oj)  in  the  form 


dds  ^ 

+  QrO^s  + 


dal. 


“1“  ^ as^rtQ  “1“  ^as^s 


Oil 


(4.6) 
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where  gniu;),  ras(a;)  are  the  Raman  gain,  EIT  profile  coefficients  respectively,  Ks(u;) 
and  Kas((^)  are  Stokes  and  anti-Stokes  mode  coupling  coefficients.  In  Eq.  4.6  fa.(z,u;) 
are  the  renormalized  Langevin  noise  operators:  fai{z,uj)  =  \/ Njc  x  Fai{z,uj).  The 
sum  is  taken  over  the  relevant  Langevin  noise  operators  {/21,  /24,  /31,  /34}- 
The  general  solution  of  Eq.  (4.6)  can  be  written  as 


(  aAL)\ 

VL(L)) 
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(4.7) 
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where  M  = 

For  the  following  derivation,  let  us  define  the  coefficients  of  matrix 


as 


Bi\  ^  -ml 
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(4.8) 


Due  to  the  linearity  of  Eq.  (4.7),  the  unknown  variables  al{L)  and  aas(O)  of  the 
backward  wave  problem  can  be  written  as  a  linear  combination  of  the  initial  boundary 
values  and  the  noise  terms: 


where 
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(4.9) 


'  A  -  -BiCi 
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The  coefficients  A,B,C  and  D  are  the  functions  of  a;,  whereas  PaiiQai  are  the 
functions  of  uj  and  2;. 
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4.4  Characteristics  of  generated  photons 


4.4.1  Stokes  and  Anti-Stokes  photon  generation  rates 

We  evaluate  the  output  generation  rates  of  Stokes  and  anti-Stokes  photons  into  a 
single  transverse  mode,  for  example,  into  a  pair  of  mode-matched  optical  hbers.  The 
generation  rates  at  the  corresponding  boundaries  zj  are 

where  subscript  j  denotes  either  Stokes  or  anti-Stokes  photon,  Zg  =  L  and  Zas  =  0. 

Of  importance  are  the  spectral  properties  of  the  generated  photons.  The  power 
spectrum  of  the  output  Stokes  and  anti-Stokes  fields  are  related  to  their  first  order 
coherence  functions  =  (aj(zj,t)aj(zj,t  +  t))  as 

H-oo 

Rj(u;)  =  ^  J  dTP^^Gf\T)  (4.13) 

—  CXD 

We  use  the  solutions  for  as{L,  uj)  and  al^(0,  — cn)  held  operators,  given  by  Eq.  (4.9), 
the  commutation  relations  for  the  input  held  operators  [aj{zj,ijj),  a^j{zj,  —oj')]  =  Lj (27rc)5(a;-l- 
o;')  and  Eq.  (4.41).  We  apply  inverse  Fourier  transformation  a{t)  =  J  duje~^‘^^a{uj) 
and  a^(t)  =  f  duje~^‘^^a^  {—uj)  and  obtain  the  Stokes  and  anti-Stokes  generation  rates 
from  Eq.  (4.12)  in  the  form 
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(4.14) 


(4.15) 


The  integrands  of  Eq.  (4.14)  and  Eq.  (4.15)  are  Stokes  and  anti-Stokes  spectral 
generation  rates  respectively. 

As  seen  from  Eq.  (4.14)  and  Eq.  (4.15),  Stokes  and  anti-Stokes  spectral  generation 
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rates  consist  of  parametric  counts,  characterized  by  the  transfer  functions  ^(a;)  and 
B(u),  and  of  noise  counts,  that  originate  from  Langevin  noise  fluctuations.  If  the 
pump  is  weak  and  far  detuned,  B(uj)  =  —C(ll’),  and  therefore  the  parametric  terms 
in  Stokes  and  anti-Stokes  rates  are  equal.  We  also  note  that  the  contribution  of 
Langevin  noise  fluctuations  to  spectral  generation  rates  can  not  be  neglected  and 
becomes  dominant  at  low  optical  depth  of  the  atomic  sample. 

Using  the  commutator  conservation  relation  for  the  Stokes  field  at  the  right  bound¬ 
ary  z  =  L,  obtained  from  Eq.  (4.9),  the  Stokes  spectral  generation  rate  can  be  ex¬ 
pressed  as 


L 

RM  =  \A\‘  -  1  +  E  / (4-16) 

g 

When  the  pump  is  far  detuned  from  the  atomic  transition  and  Stokes  photon  losses 
are  small,  the  contribution  of  Langevin  noise  fluctuation  in  Eq.  (4.16)  is  negligible. 
Thus,  the  Stokes  spectral  generation  rate  can  be  written  as 

RM  ~  |41|'  -  1  (4.17) 

Eq.  (4.17)  can  be  interpreted  in  terms  of  the  quantum  theory  of  linear  amplifica¬ 
tion  [47] :  in  the  absence  of  the  Stokes  input  beam  and  losses  for  the  Stokes  photon,  the 
Stokes  generation  rate  is  just  the  additive  noise  caused  by  the  amplification  process, 
that  is  characterized  by  the  gain  coefficient  A{u!). 

4.4.2  Two-photon  intensity  correlation  function  and  bipho¬ 
ton  wave- function 

In  order  to  address  another  important  issue  -  the  time  correlation  properties  of  the 
generated  photons,  we  calculate  the  Glauber  two-photon  correlation  function  of  time 
delay  r  between  Stokes  and  anti-Stokes  photons: 


GAsM  =  +  T)a^,(0,t  +  T)a,(L,t)) 


(4.18) 
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Using  the  Stokes  and  anti-Stokes  operators  at  the  boundaries  derived  earlier 
[Eq.  (4.9)]  and  the  commutator  relations  for  the  input  helds,  we  obtain  the  inten¬ 
sity  correlation  function  as 


Gfi,(T)  = 


X  (dl(L,-uJi)dl^(0,-U2)das(0,uj3)as(L,U4))  (4.19) 


The  intensity  correlation  function  contains  the  fourth-order  Langevin  noise  cor¬ 
relations.  According  to  Gaussian  moment  theorem  [48,  49]  they  can  be  decomposed 
preserving  the  order  into  the  sum  of  the  products  of  second  order  Langevin  noise 
correlations.  As  a  result  the  intensity  correlation  function  can  be  simplihed  to 


GAsM  =  (A(0,t  +  T)a^,(0,t  +  T))(al(L,t)a,(L,t))  +  |(a„(0,  i  +  tK(L,  «))|’ 

=  GP(0)  X  G«(0)  +  (4.20) 

The  hrst  term  in  Eq.  (4.20)  represents  flat  uncorrelated  background,  the  second 
term,  expressed  through  <l>5_as(r)  function,  describes  the  correlation  part.  <hs_a5(r) 
is  equal  to 


<l>5_as(r)  =  (das(0,t  +  T)ds(L,t)} 

L 


L 

2710 


duj  e*‘ 


bd'  +  YJ  ‘‘-QiP.i.Po 


Qi,aj  g 


(4.21) 


We  note  that  e  represents  a  two  photon  wavefunction  or  a 

(2) 

biphoton  wavefunction  on  the  condition  that  peak  value  of  normalized  g^das  » 
where  gflasi^)  =  G^lasir) / {Gf\C)  x  ^1*^(0)). 

Similarly,  we  can  dehne  <has_s(r)  =  (as(L, f)aas(0,  f  -|-  r)).  This  function  can  be 
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obtained  in  the  form 


(4.22) 


L 


27rc 


L 


For  a  wide  range  of  input  parameters  we  numerically  verify  that  <has-s(r)  = 


s—as 


For  the  case  where  Qp/Aun  «  1  and  Aa;i4/7i4  >>  1,  we  numerically  verify  that 
the  contribution  of  Langevin  noise  fluctuations  to  $as-s('r)  is  negligible,  therefore 
Eq.  (4.22)  can  be  simplihed  to 


(4.23) 


4.5  Ideal  Spontaneous  Parametric  Down  Converter 

Before  we  proceed  to  the  discussion  of  the  interesting  cases  of  the  atomic  correlation 
functions  and  photon  spectral  densities,  we  want  to  make  an  analogy  to  the  well- 
known  parametric  down  converter  in  crystals  [40].  We  consider  the  ideal  model  - 
non-degenerate  parametric  down  converter  in  which  a  generated  signal  photon  has  a 
very  slow  group  velocity  Vg  as  compared  to  an  idler  photon.  We  assume  that  both  idler 
and  signal  photon  escape  SPDC  without  losses,  therefore  the  Langevin  noise  terms  in 
Eq.  (4.6)  can  be  neglected.  In  crystals  the  coupling  coefficient  can  be  approximated  as 
a  constant  over  the  broad  spectral  range  ^^(a;)  =  Kas{uj)  =  k.  We  also  neglect  Raman 
gain  gnioj)  and  approximate  EIT  profile  as  F(a;)  =  —iuj/Vg.  Under  these  assumptions 
the  signal  and  idler  photons  have  identical  spectral  characteristics  and  rates.  Using 
the  Eq.  (4.14),  we  obtain  the  photon  spectral  density  R{ijj)  and  spectrally  integrated 
generation  rate  R  =  f  duR(u)  in  the  form: 


CHAPTER  4.  THEORY  OF  EIT  BASED  PAIRED  PHOTON  GENERATION  40 


R(uj)  =  kpL^sinc^  (  — —  ) 

\‘^VgJ  (424) 

R  =  Vg\K\‘^L. 

As  seen  from  Eq.  (4.24),  the  spectral  bandwidth  of  the  SPDC  in  crystals  is  limited 
by  Ao;  ~  27rVg/L  due  to  the  phase  mismatch  of  the  off-centered  counts  resulting  from 
the  group  delay  in  the  signal  channel.  In  principle,  the  bandwidth  can  be  made 
very  small  by  making  the  group  delay  Tg  =  L/Vg  large.  For  the  atomic  system  k  is 
proportional  to  the  atomic  density  M  =  N/V,  therefore  the  spectral  density  of  the 
generated  photons  scales  as  (AfL)^.  Taking  into  account  that  the  EIT  induced  group 
velocity  of  a  wave  is  I4  =  fl^/ {2'ji3J\fai3),  where  (T13  =  pi3a;i3/(ceo^7i3)  is  the  atomic 
cross  section  of  the  |1)  — |3)  transition,  the  total  count  rate  scales  linearly  with  AfL. 


Figure  4.2:  Normalized  signal-idler  intensity  correlation  function  for  the  ideal  SPDC 
with  the  condition  of  a  large  group  delay  LjVg  for  a  signal  photon.  R  is  the  paired 
photon  generation  rate. 

Fig.  4.2  shows  the  normalized  signal  idler  intensity  correlation  function.  Since 
the  probability  of  emitting  a  photon  pair  is  uniformly  distributed  along  the  crystal 
of  length  L  and  the  signal  photon  has  a  group  delay  relative  to  the  idler,  the 
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is  expected  to  be  the  shifted  rectangle  with  the  width  equal  to  the  group  delay  Tg  = 
L/Vg.  By  varying  the  group  delay  we  can  control  the  width  of  the  waveform.  The 
peak  height  of  is  equal  to  1/(kL)^  and  can  be  interpreted  as  a  duty  ratio:  = 

1  / {RTg).  In  the  regime  of  very  small  parametric  gain  kL  <<  1,  which  can  be  achieved, 
for  example,  by  lowering  the  intensity  of  the  pump,  the  atomic  system  can  produce 
highly  correlated  photon  pairs  gl_^  >>  1. 

4.6  EIT  based  Paired  Photon  Generator 

Now  we  turn  from  the  discussion  of  the  ideal  SPDC  to  the  discussion  of  the  EIT 
based  paired  photon  generator.  In  the  EIT  based  generator,  coupling  between  Stokes 
and  anti-Stokes  modes  is  bandwidth  limited.  Moreover,  the  generated  Stokes  photon 
undergoes  Raman  gain  gR^uj),  whereas  a  paired  anti-Stokes  photon  propagates  slowly 
and  undergoes  absorbtion  at  the  poles  of  EIT  prohle  oj  =  ±f2c/2. 

We  hrst  obtain  the  coefficients  of  Eq.  (4.6)  with  the  assumption  that  the  pump  is 
weak  and  far  detuned  from  the  |1)  — >  |4)  transition  and  Afc  =  {kp+kc  —  ks  —  kas)-z  =  0: 
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r  = 

-L  n  s 


9r  — 


R.9  —  Rn.R  — 


2iM(T'yiz{uj  +  ^712) 

G{uj) 

\  iA/’a7i3(a;  +  *713) 

W)  ’ 

\  *A/'cr7i3 


-( 


V2Awi4/  G{uj) 


(4.25a) 

(4.25b) 
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(4.25i) 

(4.25j) 

(4.25k) 


where  M  is  an  atom  density,  G(u;)  =  |r2cP  —  4(a;  +  i7i2)(a;  +  i7i3),  a  is  the  absorption 
cross  section  for  all  allowed  transitions:  a  =  au  =  a24  =  cr23  =  (J13. 
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4.6.1  Ground  state  approximation 

The  approximation  that  the  atomic  population  remains  in  the  ground  state  an  =  1 
gives  two  signihcant  diffusion  coefficients  T>i2,2i  =  2712  and  T>i3,3i  =  r3  correspond¬ 
ing  to  -^21,  -^31  Langevin  noise  operators.  Due  to  small  atomic  population  in  the 
excited  states,  the  rest  of  the  diffusion  coefficients  are  approximated  as  zeros  and  the 
corresponding  Langevin  noise  operators  are  neglected  [36]. 

We  numerically  examine  the  emission  rates  and  intensity  correlation  function  for 
the  EIT  based  paired  photon  emitter.  We  take  other  parameters  similar  to  those 
of  a  Rb  MOT:  atom  density  Af  =  10^^  atoms  per  cm^,  atomic  cross  sections  a  = 
<713  =  <7i4  =  (724  =  10“®  cm^  and  dephasing  rates  equal  to  one  half  of  the  Einstein  A 
coefficient,  i.e.,  713  =  714  =  724  =  723  =  1.79  x  10^  radians.  We  choose  the  strength 
and  the  detuning  of  the  pump  laser  from  the  |1)  — >■  |4)  transition  as  Aujn  =  24713 
and  Dp/ Nujn  =  0.1. 
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Figure  4.3:  Transmission  and  phase  mismatch  as  functions  of  detuning  cu 

Fig.  4.3  shows  the  EIT  transmission  prohle  and  phase  mismatch  as  a  function  of 
the  detuning  of  the  anti-Stokes  frequency  ui.  We  take  Dc  =  6713,  MaL  =  11  and 
7i2  =  0.  Fig.  4.4  shows  the  prohles  of  the  coupling  coefficient  |K(a;)|  and  the  Raman 
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Figure  4.4:  Coupling  constant  and  Raman  gain  as  functions  of  detuning  a;. 

gain  coefficient  Re{gR{u}))  for  the  same  parameters  as  in  Fig.  4.3. 

EIT  based  spontaneous  emitters  allow  a  variable  emission  bandwidth.  The  band¬ 
width  and  profile  of  the  spontaneous  emission  rates  are  controlled  by  the  strength  of 
the  coupling  laser  and  the  optical  depth  of  the  atomic  sample  J\faL.  By  reducing 
the  strength  of  the  coupling  laser,  the  emission  linewidth  can  be  made  much  smaller 
than  the  natural  linewidth,  with  the  minimum  width  ultimately  limited  by  the  de¬ 
phasing  rate  712  of  the  |1)  — >  |2)  transition.  At  small  optical  depth  AfaL  <  1,  the 
Raman  gain  coefficient  gR^uj)  and  the  coupling  coefficients  determine 

the  emission  spectrum.  At  high  optical  depth  J\faL  >  1,  the  EIT  transmission  win¬ 
dow  and  the  phase  mismatch,  introduced  by  a  large  group  delay  in  the  anti-Stokes 
channel,  affect  the  spontaneous  emission  spectrum. 

In  Fig.  4.5,  Fig.  4.6  and  Fig.  4.7  we  show  the  variations  of  the  coincidence  count 
rate  Rc{t)  in  a  1  ns  bin  and  the  corresponding  Stokes  power  spectral  density  depend¬ 
ing  on  the  optical  depth.  With  a  bin  size  AT  =  1  ns  much  smaller  than  the  correlation 
time,  the  coincidence  count  rate  is  obtained  from  the  intensity  correlation  function 
as  Rc{t)  =  AT{c/ L)'^G^as-s{T)-  Compared  to  the  ideal  case,  described  earlier,  the 
intensity  correlation  function  and  the  emission  spectrum  for  the  EIT  based  atomic 
SPEC  show  some  interesting  features.  The  shape  of  the  intensity  correlation  function 
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Figure  4.5:  The  oscillatory  regime:  (a)  Coincidence  count  rate  in  a  1  ns  bin  and  (b) 
Stokes  spectral  generation  rate.  J\faL  =  0.3,  hlc  =  6713  and  712  =  0. 
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Figure  4.6:  The  group  delay  regime:  (a)  Coincidence  count  rate  in  a  1  ns  bin  and  (b) 
Stokes  spectral  generation  rate.  J\faL  =  20,  hlc  =  6713  and  712  =  0. 
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Figure  4.7:  The  group  delay  regime:  (a)  Coincidence  count  rate  in  a  1  ns  bin  and  (b) 
Stokes  spectral  generation  rate.  J\faL  =  200,  flc  =  6713,  Dp  =  0.8713  and  712  =  0. 
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Figure  4.8:  Intensity  correlation  function  with  the  correlation  time  smaller  than  the 
spontaneous  decay  time,  Tj.  <  Tg  <  1/(2713).  J\faL  =  200,  lie  =  65713,  712  =  0  and 
1/(2713)  =  25  ns. 
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and  the  emission  profile  depend  on  the  relation  of  three  characteristic  times.  The  hrst 
is  the  inverse  Rabi  frequency  of  the  coupling  laser  =  27r/ the  second  is 
the  group  delay  between  Stokes  and  anti-Stokes  photons  Tg  =  L/Vg  =  2'^i‘iN'aL/D‘l 
and  the  third  is  a  minimal  pulse  length  required  to  pass  through  EIT  medium  [50] 
ln(2)7i3A/A/'crL/hlc. 

When  the  EIT  effect  is  small,  which  occurs,  for  example,  at  low  optical  depth,  the 
atomic  system  behaves  like  a  single  atom  [51].  In  such  a  regime  the  intensity  corre¬ 
lation  function  reveals  the  damped  Rabi  oscillations  (Fig.  4.5).  The  oscillations  ob¬ 
served  in  the  intensity  correlation  function  have  the  time  period  of  =  2tt/ —  7^3 
and  occur  on  the  condition  that  Tr  >  Tg,  Tp(^rnin)  and  the  coupling  laser  is  strong  enough 
to  force  the  oscillations  to  overcome  damping  flc  >  7i3-  Once  the  metastable  state  |2) 
is  excited  by  the  Raman  process  |1)  — >•  |4)  — >•  |2),  the  probability  amplitude  between 
1 2)  and  |3)  oscillates  due  to  the  strong  interaction  of  the  atoms  with  the  resonant 
coupling  beam. 

In  Fig.  4.6  we  show  the  intensity  correlation  function  and  Stokes  emission  spectrum 
in  the  group  delay  regime,  where  Tg  >  Tr  and  Tg  >  Tp^rnin)-  In  this  regime  the  width 
of  the  intensity  correlation  function  is  approximated  by  Tg.  Moreover,  the  frequency 
range  over  which  the  spontaneous  generation  occurs  is  hltered  by  the  EIT  window 
and  mostly  controlled  by  the  phase-matching  in  the  presence  of  large  group  delay 
in  the  anti-Stokes  channel.  A  sufficiently  wide  EIT  window  Tg  >  Tp^rnin)  requires 
high  optical  depth  AfaL  >  10.  In  the  presence  of  the  non-zero  dephasing  712,  the 
maximum  group  delay  and  therefore  the  maximum  width  of  the  correlation  function 
is  limited  by  1/712- 

By  increasing  the  optical  depth  of  the  atomic  sample,  the  EIT  window  can  be  made 
substantially  larger  than  the  emission  bandwidth.  A  large  EIT  window  might  be  very 
useful  for  such  applications  as  bi-photon  waveform  control  and  shaping.  Fig.  4.7 
shows  the  intensity  correlation  function  and  the  Stokes  emission  spectrum  under  such 
a  condition,  where  AfaL  =  200.  The  tail  of  the  correlation  function  decays  on  the 
time  scale  of  Tp(^rnin)- 

We  also  note  that  at  high  optical  depth  it  becomes  possible  to  achieve  the  Stokes- 
anti-Stokes  correlation  function  with  the  width  shorter  than  the  spontaneous  decay 
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time  1/(2713).  Fig.  4.8  shows  such  a  correlation  function  that  is  obtained  for  NaL  = 
200  and  =  65713,  where  Tr  <  Tg  <  1/(2713). 

We  turn  next  to  the  sharp  peak  at  the  leading  edge  of  the  correlation  function  of 
Fig.  4.7.  This  sharp  peak  is  Sommerfeld-Brillouin  precursor  [52].  It  has  been  observed 
in  our  theoretical  plots  for  a  long  time.  But  only  recently,  when  it  was  experimentally 
measured,  we  got  interested  in  its  nature.  We  thank  Daniel  J.  Gauthier  who  suggested 
that  this  peak  is  the  result  of  simultaneously  generated  Stokes  and  anti-Stokes  photons 
that  travel  at  nearly  the  speed  of  light  in  vacuum  and  arrive  near-simultaneously  at 
the  photodetectors.  He  also  pointed  out  to  the  similarity  of  this  peak  to  precursors 
which  have  been  extensively  studied  [52,  53,  54].  A  precursor  has  an  approximate 
width  that  is  equal  to  the  opacity  width  of  the  atomic  transition  in  the  optically  thick 
medium.  Similar  to  precursors,  the  sharp  leading  edge  peak  has  a  opacity  width  of 
the  EIT  profile  [55]. 

4.6.2  No  ground  state  approximation 

If  the  pump  is  weak  and  far  detuned  and  therefore  most  of  the  atomic  population  is 
in  the  ground  state,  we  verify  that  the  ground  state  approximation  gives  a  correct 
prediction  for  the  Bi-photon  function  and  the  Stokes  generation  rate.  Nevertheless 
it  does  not  properly  account  for  an  atom  return  to  the  ground  state  |1).  Ideally  one 
would  expect  the  Stokes  and  anti-Stokes  rates  to  be  equal,  since  an  atom,  making  a 
complete  cycle  on  the  energy  level  diagram  (Fig.  5.7),  returns  to  the  ground  state. 
Even  at  zero  dephasing  rate  712  =  0  of  the  |1)  — >•  |2)  transition,  the  ground  state 
approximation  predicts  the  anti-Stokes  generation  rate  to  be  smaller  than  the  Stokes 
generation  rate.  Eor  example,  Ras/Rs  =  0.65  for  AfaL  =  10,  f2p/Au;i4  =  0.1  and 

=  6713. 

In  order  to  treat  properly  an  atom’s  return  to  the  ground  state  we  will  retain  in 
Eq.  (4.6)  all  four  Langevin  noise  operators  {/21, /24, /si, /34}  and  take  into  account 
small  incoherent  population  in  excited  states,  resulting  from  the  steady  state  solutions 
of  Eqs.  (4.27a)-(4.27f).  With  these  inclusions  the  solution  of  Eq.  (4.6)  predicts  the 
Stokes  and  anti-Stokes  spectral  generation  rates  to  be  equal  at  a  zero  dephasing  rate 
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7i2  =  0  [Fig.  4.9(a)].  A  non-zero  depliasing  rate  712  7^  0  reduces  EIT  and  therefore 
introduces  additional  losses  for  anti-Stokes  photons.  As  a  result  the  output  Stokes 
rate  exceeds  the  anti-Stokes  rate.  For  example,  for  AfaL  =  10,  Dp/Auji^  =  0.1, 
=  6713  and  7i2  =  O.6713  we  obtain  Ras/Rs  ~  0.8.  The  corresponding  Stokes  and 
anti-Stokes  spectral  generation  rates  are  shown  in  Fig.  4.9(b).  We  believe  that  the 
“missing”  anti-Stokes  photons  are  absorbed  and  are  then  reemitted  in  a  solid  angle 
of  dvr.  The  atomic  sample  is  assumed  to  be  optically  thin  in  the  radial  direction. 


Figure  4.9:  Stokes  (solid  curve)  and  anti-Stokes  (dash  curve)  spectral  generation  rates 
at  (a)  zero  dephasing  712  =  0  and  (b)  non  zero  dephasing  712  =  O.6713.  AfaL  =  10, 
flc  =  6713  flp  =  2.4713,  Aa;i4  =  24713 

Compared  to  the  ideal  SPDC,  where  each  generated  signal  photon  has  its  paired 
idler  photon,  the  real  atomic  system  has  uncorrelated  noise  counts  in  both  Stokes 
and  anti-Stokes  channels  that  result  from  Langevin  noise  fluctuations.  In  Fig.  4.10 
we  examine  the  dependence  of  the  Stokes  (anti-Stokes)  and  paired  count  rates  on 
the  optical  depth  AfaL.  The  parameters  for  the  curves  are  712  =  0,  flc  =  5713, 
Aa;i4  =  24713  and  fIp/Aa;i4  =  0.1.  The  paired  count  rate  {Rp)  is  dehned  as  the  area 
under  the  Stokes-anti-Stokes  coincidence  count  rate  function  minus  the  area  under 
the  uncorrelated  background.  One  may  show  that  Rp  ~  l/(27r)  J  duj  At  small 

optical  depth  the  paired  rate  scales  quadratically  with  the  optical  depth  and  is  much 
smaller  than  the  Stokes  rate.  At  high  optical  depth  the  paired  rate  varies  linearly 
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Optical  Depth 


Figure  4.10:  Stokes,  anti-Stokes  (solid  curve)  and  Paired  (dashed  curve)  photon  gen¬ 
eration  rates  as  a  function  of  the  optical  depth,  flc  =  5713,  712  =  0,  Aa;i4  =  24713  and 
flp/Acuii  =  0.1.  At  optical  depth  of  100,  the  paired  rate  reaches  90%  of  the  Stokes 
rate. 

with  MaL  and  converges  logarithmically  to  the  Stokes  emission  rate. 

In  this  chapter  we  describe  the  theory  of  paired  photon  generation  in  double- 
A  atomic  system.  With  low  parametric  gain  and  high  optical  depth  we  show  that 
the  system  can  produce  highly  correlated  photon  pairs.  The  shape  of  the  intensity 
correlation  function  and  the  emission  bandwidth  depend  on  the  coupling  laser  Rabi 
frequency  and  the  optical  depth  of  the  atomic  sample.  Compared  to  the  ideal  SPDC, 
paired  photon  generation  in  the  double-A  atomic  system  is  affected  by  Raman  gain  in 
the  Stokes  channel  and  EIT  in  the  anti-Stokes  channel.  EIT,  through  the  absorption 
at  the  poles,  cuts  the  emission  bandwidth.  In  order  to  enter  a  regime  where  the  EIT 
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window  is  sufficiently  large  and  therefore  the  emission  bandwidth  is  controlled  to  a 
large  extent  by  the  phase-matching  process  in  the  presence  of  large  group  delay,  the 


optical  depth  of  the  atomic  sample  has  to  be  large  AfaL  >  10.  High  optical  depth 


substantially  reduces  the  influence  of  Langevin  noise  fluctuations  on  paired  photon 
generation  so  that  the  Stokes  and  anti-Stokes  photons  are  generated  mostly  in  pairs. 
We  therefore  suggest  the  use  of  a  cigar  shaped  atomic  cloud  with  high  optical  depth 
in  the  longitudinal  direction. 

4.7  Mathematical  Apparatus 

4.7.1  Collective  slowly  varying  atomic  operators 

To  describe  the  quantum  properties  of  the  atomic  system  we  use  the  collective  slowly 
varying  atomic  operators  [56,  57,  58,  46,  59]  ajk^zR),  defined  as 


(4.26) 


where  the  averaging  is  done  over  each  atom  i  in  a  small  interval  Az  that  contains 
large  number  of  atoms  W  >>  1.  The  slowly  varying  variables  are  assumed  to  stay 
unchanged  over  A^;.  z/41  =  —  oji  +  Aa;i4,  1/42  =  a;4  —  a;2  -|-  Aa;i4,  1/44  =  —1^41, 

z/24  =  -z/42,  the  rest  of  Vjk  =  - ^k-  hi  =  hs  ■  z,  ^42  =  ks-z,  kn  =  kp-z,  k^2  =  kc-z 

are  the  projections  of  the  anti-Stokes,  Stokes,  pump  and  coupling  /c-vectors  on  the  2; 
axis,  /C43  =  kii  —  /C31,  ^21  =  —^24  +  kii,  the  rest  of  kjk  =  —kkj. 


4.7.2  Heisenberg-Langevin  equations 

The  full  set  of  Heisenberg-Langevin  equations  for  the  four  state  system  consists  of  16 
equations.  Here,  we  show  explicitly  10  of  them,  the  other  6  for  the  adjoint  off-diagonal 
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atomic  operators  are  not  shown. 


7w0'44  —  -^44  —  r4Cr44 

Ot 

+l  I  “^<^41  +  5'sOsCr42  —  (7sa^cr24 - ^*^14 

d  _  ~ 

^CT4i  =  F41  —  (741  +  iAu;i4)  (741 

at 

f  D*  ^ 

Ti  (  — 5'ash]js(743  +  gs%<y 21  +  “^  (<^11  ~  <^44) 
^CTss  =  F33  —  r3a‘33 

f  ^  Vl* 

Ti  igas(IasO'31  +  “^^32  “  S'ashLcnS - ^^^3 

—  (732  =  F32  —  7320'32 

ot 

(  __  D* 

+i  f  -gasaisai2  +  gsala^i  +  (^33  -  ^22) 

—  a'22  =  F22  +  r32a'33  +  r42a44 

I - ^*^32  ~  5'sflsO'42  +  fi'sfls'^24  +  “^*^23 

d  _  ~ 

—  (7ii  =  Fii  +  r3i(733  +  r4i(744 

ot 

D 

-gas^asd'31 - 5'as®L^13  +  “^^14 


(4.27a) 


(4.27b) 


(4.27c) 


(4.27d) 


(4.27e) 


(4.27f) 


+i 
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^(T34  =  F34  —  (734  —  l/XtOip  (J34 

[Tt  __  __  D* 

pi  (  “^*^31  +  5'sfls5‘32  —  - ^^24 

^  ~  -T 

— (T31  —  i^3i  —  7310-31 

at 

.  /  .t  /-  ~  ^ 

+1  I  fi'asOL  (<^33  ~  <^ll) - ^*^21  +  “^*^34 

d  _  ~  . .  X  _ 

^0-24  =  -^24  —  (724  —  *Au;i4)  0^24 

+i  ^”^^21 - ^^34  +  Qs^s  {^22  —  0-44) 
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Ot 

+*  I - ^*^31  ~  5's®s<^41  +  gasCas^23  +  “^^24 


(4.28a) 


(4.28b) 


(4.28c) 


(4.28d) 


Here,  for  simplicity,  we  assume  that  Ak  =  {kp  Pkc  —  ks  —  kas)  -2  =  0.  In  Eq.  (4.27a)- 
Eq.  (4.28d)  E,  is  the  total  decay  rate  from  state  |z),  Tij  is  the  decay  rate  from  state 
|i)  to  state  \j)  and  ■jij  is  the  dephasing  rate  bet-ween  state  |i)  and  state  \j).  The 
dephasing  rates  for  the  double-A  system  in  the  absence  of  the  collisional  dephasing 
can  be  obtained  from  total  decay  rates  r3  and  Fi  from  state  |3)  and  |4)  to  t-wo  ground 
states  |1)  and  |2)  as 


r4 

741  —  742  — 


r3  -i-  r4 


743  = 


2 


(4.29) 

(4.30) 

(4.31) 
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4.7.3  Linearization  procedure 

In  zeroth  order  perturbation  expansion,  in  which  and  Gas  go  to  zero,  the  Heisenberg- 
Langevin  equations  for  dn,  5'i4,  d22,  5'23,  5'32,  dsa,  a'4i,  a'44  atomic  operators  are  decou¬ 
pled.  Under  the  assumption  that  pump  and  coupling  beams  propagate  without  de¬ 
pletion,  we  obtain  the  steady  state  solution  for  Eqs.  (4.27a)-(4.27f)  in  form 

^jk  =  (^jk)  +  XI  (4.32) 

With  the  dehnition  of  the  denominator  as 

T  =  r3i  {Tl  +  AAul  +  2\Dp\^)  Iff/ 

+  r42(r2  +  2|ff,|')  (4.33) 

the  steady  state  expectation  values  for  the  zeroth  order  atomic  operators  are  equal 
to 


/-o  \  ^31  ^  (r4  -|-  4Aa;^4  -|-  hip 

Wii/  “  2" 

(4.34a) 

r42  (rl  +  p/) 

W22/  -  rjn 

(4.34b) 

/~0  \  ^42  hichip  ^ 

W33/  -  j, 

(4.34c) 

/~0  \  ^31  hichip  ^ 

\Cr44)  ^ 

(4.34d) 

/_o  \  Esi  (2Aci;i4  —  ir4)  hic  ^ hip 

(‘^14/  rj^ 

(4.34e) 

/_o  \  *r3r42hic  hip  ^ 

('^23/  ~  rp  ) 

(4.34f) 

In  the  hrst  order  expansion,  we  substitute  the  zeroth  order  solution  for  the  atomic 
operators  Eq.  (4.32)  into  the  remaining  Heisenberg-Langevin  equations  for  d2i,5'24, 
d3i,  d34  and  their  adjoint.  Neglecting  higher  order  terms  like  emnFmn^s  or  emnFmn^\s 
we  obtain  the  linearized  equations.  We  note  that  the  linearized  Eqs.(4.28a)-(4.28d) 
for  a‘21,  d24,  d3i,  d34  represent  an  independent  set  of  equations  and  can  be  decoupled. 
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For  clarity  we  write  them  in  vector  form 

— ai  =  +  A^a  +  Fi  (4.35) 

where  ai  =  {0-21, 0-24,  osi,  0-34},  Fi  =  {E21,  E24,  H31,  F34},  a  =  {ds,dlJ,  matrix  A 
depends  on  the  dephasing  rates  7jfc  and  pump  and  coupling  laser  rabi  frequencies 
Qp,  Dc,  matrix  M  depends  on  the  zeroth  order  solution  for  the  atomic  operators 

{^%)- 

4.7.4  Langevin  noise  operators  and  their  diffusion  coefficients 

By  analogy  with  the  collective  slowly  varying  atomic  operators  djk{z,  t),  the  collective 
Langevin  noise  operators  are  dehned  as 

=  ^  E  (4.36) 

We  assume  that  a  Langevin  noise  operator  for  a  single  atom  is  5-correlated  so  that 

T2mDA))  =  (4.37) 

where  (...)  denotes  the  average  over  the  reservoir,  atomic  diffusion 

coefficient  for  an  atom. 

Now  we  consider  the  second  order  correlations  for  the  collective  Langevin  noise 
operators 

{F„{t,z)F,,e{t',z'))  =  A  E  (4.38) 

Introducing  the  average  atomic  diffusion  coefficient 

FjkJ'k'it-,  z)  =  '^jkj'k'i^) 

^  ZiGN, 


(4.39) 
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the  noise  correlations  can  be  expressed  as 

-  t')d(z  -  z!)  (4.40) 

In  case  Djkj'k'  {t,  z)  is  independent  of  t,  in  frequency  domain  the  noise  correlations 
are 

{Fjk{ijJ,z)Ffk'{i^\z'))  =  ——Vjk,fkO{LJ  +  uj')6{z  -  z')  (4.41) 

The  diffusion  coefficients  Djkj'k'  can  be  obtained  from  the  Heisenberg-Langevin 
equations  (4.27a)-(4.28d)  using  the  generalized  fluctuation-dissipation  theorem  [60, 
61].  Here  we  show  the  diffusion  coefficients  for  the  Langevin  noise  operators  of  interest 
Ai,  -^24,  hi,  -^34  and  their  adjoint  F12,  F42,  hs,  hs 
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0 
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(5-44}  Ta  y 

(4.43) 


where  a*  denotes  {21, 24, 31, 34}  subspace  for  the  atomic  operators,  denotes 


(12,42, 13,43}  subspace  for  the  adjoint  atomic  operators. 


